D-wave correlated Critical Bose Liquids in two dimensions 
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We develop a description of a new quantum liquid phase of interacting bosons confined in two 
dimensions which possesses relative D-wave two-body correlations. We refer to this stable quantum 
phase as a D-wave Bose Liquid (DBL). The DBL has no broken symmetries, supports gapless boson 
excitations which reside on "Bose surfaces" in momentum space, and exhibits power law correlation 
functions characterized by a manifold of continuously variable exponents. While the DBL can be 
constructed for bosons moving in the 2d continuum, the state only respects the point group sym- 
metries of the square lattice. On the square lattice the DBL respects all symmetries and does not 
require a particular lattice filling. But lattice effects do allow for the possibility of a second distinct 
phase, a quasi-local variant which we refer to as a D-wave Local Bose Liquid (DLBL) . Remarkably, 
the DLBL has short-range boson correlations and hence no Bose surfaces, despite sharing gapless 
excitations and other critical signatures with the DBL. Moreover, both phases are metals with a 
resistance that vanishes as a power of the temperature. We establish these results by constructing 
a class of many-particle wavefunctions for the DBL, which are time reversal invariant analogs of 
Laughlin's quantum Hall wavefunction for bosons in a half-filled Landau level. A gauge theory 
formulation leads to a simple mean field theory and a suitable N— flavor generalization enables in- 
corporation of gauge field fluctuations to deduce the properties of the DBL/DLBL in a controlled 
and systematic fashion. Various equal time correlation functions thereby obtained are in qualitative 
accord with the properties inferred from the variational wavefunctions. We also identify a promis- 
ing microscopic Hamiltonian which might manifest the DBL or DLBL, and perform a variational 
energetics study comparing to other competing phases including the superfluid. We suggest how the 
D-wave Bose Liquid wavefunction can be suitably generalized to describe an itinerant non-Fermi 
liquid phase of electrons on the square lattice with a no double occupancy constraint, a D-wave 
metal phase. 
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I. INTRODUCTION 

The principal roadblock impeding progress in disen- 
tangling the physics of the cuprate superconductors is 
arguably our inability to access quantum ground states 
of two-dimensional itinerant electrons which are qualita- 
tively distinct from a Landau Fermi liquid. Overcoming 
this obstruction is of paramount importance, indispens- 
able in explaining the strange metal behavior observed 
near optimal doping and a likely requisite to account 
for the emergent pseudo-gap at lower energies. A puta- 
tive underlying paramagnetic Mott insulator provides the 
scaffolding for a popular class of theories, which view the 
pseudo-gap as a lightly doped spin liquid ,ii2i2idiii Signif- 
icant progress has been made in developing the ground- 
work and there now exists a well established theoretical 
framework to describe a myriad of distinct spin liquids. 
For the cuprates the most promising spin liquids are de- 
scribed in terms of fcrmionic spinons minimally coupled 
to a compact U(l) gauge field and moving in various 
background fluxes. Upon doping formidable challenges 
arise. Bosonic holons carrying the electron charge be- 
come mobile carriers and lead to electrical conduction. 
But at low temperatures Bose condensation appears in- 
evitable and this leads to Fermi liquid behavior, either a 
metal with conventional Landau quasiparticlcs or a BCS 
superconductor if the spinons are paired. Accessing a 



pseudo-gap or a strange metal which conduct electricity 
despite the absence of long-lived Landau quasiparticles 
requires doped holons that form an uncondensed quan- 
tum Bose fluid rather than a condensed superfluid. But 
is this possible, even in principle? If possible, what prop- 
erties would such a putative "Bose metal" exhibit?— iLSbL 
And what theoretical framework is appropriate? The 
slave-particle gauge theory approach has been stymied by 
this stumbling block for over 15 years. In this paper we 
provide (some) answers to these questions by construct- 
ing explicit examples of such unusual phases of bosons 
that may offer some routes out of the conundrum. 

Our goal, then, is to access and explore uncondensed 
quantum phases of 2d bosons which are conducting flu- 
ids but not superfluids. Specifically we have in mind 
hard core bosons moving on a 2d square lattice, but seek 
to construct states which do not require particular com- 
mensurate densities. While our construction can be im- 
plemented for bosons moving in the 2d continuum Eu- 
clidean plane, the states will only possess the reduced 
point group symmetry of the square lattice. Despite 
our interest in time reversal invariant quantum ground 
states, our technical approach will be strongly informed 
by theories of the fractional quantum Hall effect (FQHE). 
In some regards, the quantum phases that we construct 
are time reversal invariant analogs of the Laughlin state 
for bosons in a half-filled Landau level. But the physi- 
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cal properties of the phases will be dramatically different 
from the incompressible FQHE states, and will have gap- 
less excitations and metallic transport for example. 

To motivate and illustrate our approach, it will be 
helpful to briefly revisit the bosonic FQHE. Consider the 
Laughlin wavefunction for bosons in a half-filled Landau 
level, 
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Much of the physics of the Laughlin state has its ori- 
gin in the structure of zeroes, which reveals that any 
two particles upon close approach in real space are in 
a relative two-body d + id state, (Zi — Zj) 2 . Our first 
objective is to construct a time reversal invariant (real) 
wavefunction in which particle pairs arc similarly in a 
relative d xy state. A clue is offered by noting that the 
Laughlin wavefunction is the square of a Vandermondc 
determinant, \f , J ,=i/2 = (detV) 2 . In the Vandermonde 
determinant, all pairs of particles are in a relative p + ip 
state, (zi — Zj). Upon squaring, the two p + ip states 
combine to form a single d + id state, which is essentially 
just addition of angular momentum. 

This suggests constructing a time reversal invariant bo- 
son wavefunction in zero magnetic field by simply squar- 
ing a determinant constructed from momentum states 
within a Fermi sea, 

*(r 1 ,r 2 ,...,r w ) = (dete 4k '' r 2 , (S-type). (2) 

Let's examine the nodal structure. The generic behav- 
ior of each fermion determinant when any two particles 
are taken close together is dictated by Fermi statistics 
and reality of the wavefunction and has the functional 
form (r.i — Tj) ■ I. The unit vector I will depend in a 
complicated way on the location of all the other particles 
(see Ref. [13 for a discussion and illustrations of the free 
fermion nodes). When I = y, this is a, p x form vanishing 
along a nodal line (in the relative coordinate) parallel to 
the x-axis. Unfortunately, in contrast with the Laugh- 
lin case, squaring the determinant leads here to an "ex- 
tended" s-wave form with a quadratic nodal line rather 
than the desired d-wave. 

But consider instead multiplying together two fermion 
determinants, each constructed by filling up a Fermi sea 
of momentum states, but with different Fermi surfaces. 
A wavefunction with d xy two-particle correlations can be 
constructed by choosing two elliptical Fermi surfaces, one 
with its long axis along the x-axis and the other rotated 
by 90 degrees, as illustrated in Fig. [1] 

*(ri, r 2 , ...,r N ) = (det) x x (det). y , (D^-type) , (3) 

where the short-hands (det)^ and (det)j, represent the 
corresponding Slater determinants. In the limit of ex- 
treme eccentricity of the ellipses, the wavefunction will 
have the desired d xy form when two particles are brought 
close together, (xi — Xj)(yi — Uj)- Away from this limit, 



FIG. 1: DBL wavefunction Eq. (|3} is obtained by multiplying 
two Slater determinants corresponding to distinct but sym- 
metry related Fermi surfaces, which can be viewed as Fermi 
surfaces of di and (fe slave fermions. The left panel shows an 
example with closed elliptical Fermi surfaces, while the right 
panel is a case with open Fermi surfaces in the first Brillouin 
zone of the square lattice. 




FIG. 2: "Nodal picture" of the continuum DBL wavefunction 
as probed by moving a single particle (whose initial location 
is marked by the filled circle) while keeping the rest of the 
particles fixed (open circles). The signs of the wavefunction 

are indicated in several places to bring out the H 1 — pattern 

upon encircling a target particle. 



the two nodal lines will not align precisely with the x and 
y axes, but upon taking one particle around the other the 
wavefunction will exhibit the same sign structure as a d xy 

orbital, H 1 — , changing sign twice. A picture of such 

sign structure as seen by a test particle is shown in Fig. [5] 
Our principal thesis is that this wavefunction captures 
qualitative features of a new quantum liquid phase of 
bosons, which we call "D- wave-correlated Bose Liquid" 
(DBL) . But a variational wavefunction does not provide a 
complete characterization of a quantum phase, and can- 
not be used to address its stability. As for the FQHE, 
a field theoretic approach, such as Chcrn-Simons gauge 
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theory, is both desirable and ultimately necessary. Since 
we require time reversal invariance, Chern-Simons is in- 
appropriate, but a tractable field theoretic framework for 
the DBL will nevertheless turn out to be a gauge theory. 

Indeed, we follow closely the gauge theory approaches 
to spin liquids in quantum antiferromagnets,-^^ which 
after all are lattice bosonic systems. But there are some 
notable and important differences for the DBL which we 
discuss below. The variational wavefunction in Eq. j3|), 
being a product of two fermion determinants, naturally 
suggests expressing the boson creation operator as a 
product of two fermion operators, 

&t(r)=4(r)4(r), (4) 

with the 2d position r either continuous or denoting the 
discrete sites of a square lattice. In a general case, this de- 
composition introduces a local SU(2) gauge redundancy 
well-known in the slave-fermion treatments of the spin- 
1/2 antiferromagnet (see Ref . [Til for a recent comprehen- 
sive discussion) . The S-type wavefunction corresponds to 
a so-called SU(2) liquid, while the DBL wavefunctions 
are t/(l) liquids and require a C/(l) gauge field which 
is minimally coupled to the two fermions with opposite 
gauge charges (the product d\d\ is then a gauge neutral 
composite, the physical boson). 

Within slave particle theory a mean field state with a 
fixed background gauge "magnetic" flux is chosen. This 
choice is constrained by symmetry: All the physical sym- 
metries of the boson model, even if not respected by the 
mean field ansatz, must be present in the full gauge the- 
ory once gauge transformations are allowed. The col- 
lection of these symmetry transformations in the gauge 
theory is called the Projective Symmetry Group (PSG)Ai 
For lattice boson theories many different mean field states 
with different PSG's are possible, while little is offered to 
guide which mean field state is the correct one. 

In the continuum, the symmetries of the Euclidean 
plane together with time reversal invariance are more 
restrictive and we can only offer two PSG mean field 
ansatzc in this case. One is a mean field state with zero 
gauge flux and spherical Fermi seas for both fermions; 
this is an SU (2) ansatz and corresponds to the S-type 
wavefunction considered earlier. A different U{\) ansatz 
is obtained by considering d\ and efe fermions each mov- 
ing in a uniform field but of opposite signs for the two; 
in this case, the boson wavefunction is schematically 
* = (det) x (det)*, e.g., <J = |detV| 2 = |tf„ =1/2 | to 
give a concrete example; this wavefunction is again non- 
negative and has off-diagonal quasi-long-range orderi 12 i 13 

For bosons moving in continuous space but with only 
the square lattice point group symmetry, we can offer 
only four PSG mean fields; all of them have zero flux, 
but differ in the way the 90 degree rotation and mirror 
symmetries are realized. One example has both d\ and 
d,2 Fermi surfaces individually respecting the symmetries 
but otherwise independent of each other. Two more ex- 
amples are the already introduced D xy state, Eq. ^ and 
Fig. [TJ with elliptical Fermi surfaces elongated along the 



x and y axes, and a similarly constructed Y) x i_ y i state 
with ellipses along x ± y; in both cases, the d\ and di 
Fermi surfaces are related by the 90 degree rotation, but 
the two states differ in the way the mirror symmetries 
are realized. In the fourth example, each Fermi surface 
is invariant under the rotation but not under the mir- 
rors, which instead are realized by interchanging the two 
fermions. 

Each of the four presented ansatze can be loosely called 
a D-wave Bose liquid as far as the nodal pictures like that 
in Fig. [2] are concerned; the gauge theory analysis would 
also be rather similar. From now on, we will focus on the 
D xy boson liquid. The associated mean field wavefunc- 
tion once Gutz wilier projected into the physical Hilbcrt 
space, d\d\ — d\d2, is precisely the variational wavefunc- 
tion in Eq. A full incorporation of gauge fluctuations 
about the mean field performs, in principle, this projec- 
tion exactly. But in practice the gauge theory approach 
allows incorporation of slowly varying gauge fluxes and 
is different in this respect from the wavefunction. Pro- 
vided the gauge theory is not in a confined state and no 
physical symmetries are broken, an exact description of 
the putative DBL phase is thereby obtained. Herein, we 
introduce a large N generalization of the gauge theory 
which allows a controlled and systematic treatment of 
the gauge fluctuations order by order in powers of 1/N. 
We can then perturbatively extract physical properties of 
the DBL phase. Whether the DBL phase survives down 
to N = 1 is not something that we can reliably address. 

In describing the D xy liquid phase for lattice bosons 
there is considerable freedom when choosing the pre- 
cise form of the slave fermion hopping in the mean field 
Hamiltonian. The simplest choice is near neighbor hop- 
ping with different amplitudes in the x and y directions. 
But even with this restriction, there are two different 
possible Fermi surface topologies, closed or open, as de- 
picted in Fig. [TJ The Fermi surface topology has rather 
dramatic consequences for the nature of the associated 
D-wave Bose fluid. When the Fermi surfaces are open 
and have no parallel Fermi surface tangents, the result- 
ing phase has a quasi-local character - the boson Green's 
function is found to fall off exponentially in space. This 
is a distinct quantum phase which we refer to as a D- 
Wave Local Bose liquid (DLBL). The DLBL is intrinsi- 
cally very stable against gauge fluctuations and we are 
fairly confident that when present it will be a stable quan- 
tum phase. The effect of gauge fluctuations on the DBL 
phase, however, are more subtle. Within our systematic 
large N approach we do find a regime where the DBL 
is stable, but it is less clear that this will survive down 
to N = 1. The Gutzwillcr wavefunction does appear to 
describe a DBL phase with properties that are consistent 
with those inferred from the gauge theory. This provides 
some support for the stability of the DBL phase. 

Together the gauge theory and the variational wave- 
function provide a consistent and rather complete picture 
of both the DBL and the DLBL. Here we briefly high- 
light some of the important characteristics. In the DBL, 
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the single particle boson Green's function, G&(r, r) = 
(M(r, r)6(0, 0)), decays as an oscillatory power law at 
equal times, 



G 6 (r,0) 



cos[(k Fl - kf 2 ) ■ r] 

|r| 4 "" 
cos[(k Fl + k F2 ] 



r - 3tt/2] 



(5) 



where the two wave vectors k Fl 2 depend on the observa- 
tion direction f , as docs the anomalous exponent rj > 0. 
As the direction of f is rotated in real space, the wave 
vectors kj^ 2 and kp^ ± kp 2 trace out closed momentum 
space curves, and for the DBL in the continuum these 
curves will have the topology of a circle. 

If one measures the boson Green's function and finds 
it to be of the above form, it is natural to refer to the 
kpj ± kp 2 as "Bose surfaces." Within the gauge theory 
description the origin of these singular surfaces can be 
traced to the Fermi surfaces of the constituent fermions. 
Specifically, kp 1 2 are locations on the two Fermi surfaces 
where the surface normals are along the observation di- 
rection f, and the momentum space areas enclosed by 
each of these surfaces will equal (2-7r) 2 pb- Note that the 
kp 1 2 surfaces can be uniquely reconstructed from the 
measured Bose surfaces kp^ ± kp 2 . In this sense, the 
Bose surfaces in the DBL contain information analogous 
to Luttinger's Fermi surface volume theorem in the Fermi 
liquid. The Bose surface can also be extracted directly 
from the DBL wavefunction by using variational Monte 
Carlo to compute the equal time boson Green's function, 
and can also be indirectly inferred by using the shape of 
the two Fermi surfaces input into the (det)^ and (det). y 
factors. It is both reassuring and quite remarkable that 
these two coincide. 

The spatially local Green's function, Gb(0,r), in the 
DBL falls off in time as a power law which is at least as 
slow as 1/t 2 corresponding to the local boson tunneling 
density of states Ab{E) ~ E; these are mean field results 
obtained by combining two fermions each with a finite 
density of states, but we suspect the time decay may 
actually be slower upon including gauge fluctuations. 

The behavior of the boson Green's function in the 
DLBL is quite different. The time decay at r = is 
a power law, G&(0,t) ~ r~ 2 , which is an exact result 
insensitive to lattice scale details. But the equal time 
Green's function falls off exponentially in space in the 
DLBL, Gb(r, 0) ~ exp(— r/£). Despite this, a two boson 
box correlator, 



B b (x) = (tf(0,0)tf(x,x)b(x,0)b(0,x)) 



(6) 



falls off as a (non-oscillatory) power law at large distances 
in the DLBL, B(x) ~ — a: -8 , with both the sign (nega- 
tive) and the exponent being universal and insensitive to 
lattice scale physics. Paradoxically, this seems to imply 
that a pair of bosons injected into the DLBL on oppo- 
site corners of a square box can move more readily than 
a single injected boson. But in a strongly interacting 



quantum state such a single-particle interpretation can 
be misleading - the dynamics of an injected boson is not 
that of a weakly interacting quasiparticlc. 

In both the DBL and the DLBL the density-density 
correlator -D&(r) = (p&(r)p&(0)) behaves at large dis- 
tances as, 



D b (r) 



a=l,2 



cos[2k/ 
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(7) 



where again both the wave vectors k Fl 2 and the anoma- 
lous exponent "f a depend on the observation direction r. 
Measurement of such correlation can be used to extract 
k-F 12 in the DLBL, which are not accessible from the bo- 
son Green's function here. Upon rotating the unit vector 
r in real space, k^ 2 will again trace out the Fermi sur- 
faces of the underlying fermions which for the DLBL will 
be open. Both the DBL and the DLBL are conductors 
with a resistance varying with temperature as R ~ T 4 / 3 . 
This is particularly striking for the DLBL where the bo- 
son Green's function is short-ranged. Evidently, it is not 
possible to understand the metallic transport in terms of 
the motion of weakly interacting bosons. In a sense it 
is the fcrmionic constituents which transport the charge, 
but even this is not quite right since the d\ , e?2 fermions 
do not exist as well defined quasiparticlc excitations, be- 
ing strongly scattered by the gauge fluctuations. 

The lattice gauge theory formulation can be used to 
motivate a lattice boson Hamiltonian which might plau- 
sibly have the DBL/DLBL as a ground state. The sim- 
plest such model consists of a near neighbor boson hop- 
ping term supplemented by a four-site "ring" exchange 
term: 



H — Hj + H4 , 
Hj = -J ]T {b\b r+il + h.c.) 



(8) 
(9) 



Hi = Jf 4 ^(^MU?W + h - c -) > ( 10 ) 



with J, K4 > 0. This model is fully specified by two di- 
mcnsionlcss numbers, the ratio J/K and the boson filling 
< Pb < 1. This ring Hamiltonian with J > and K4 < 
was introduced and studied by Paramekanti et al.^ 
and later studied extensively with quantum Monte Carlo 
by Sandvik et al.^ Melko et al.^ and Rousseau et alrH 
There is no sign problem in this regime, and it was pos- 
sible to access large system sizes and low temperatures. 
But for J > and K4 > 0, there is a sign problem 
since the Hamiltonian does not satisfy the Marshall sign 
conditions, and one expects the ground state wavefunc- 
tion to take both positive and negative values. For this 
ring model we have evaluated the variational energetics 
for the DBL/DLBL wavefunctions and compared these 
to the energy of a superfluid wavefunction of the usual 
Jastrow form. Within this necessarily limited energetics 
study we do find a region of parameter space with small 
J/K4 and near half filling where the DLBL wavefunction 
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has the lowest energy. In view of its quasi-local charac- 
ter, we suspect that DMRG studies could be fruitful in 
helping establish whether or not the DLBL is present in 
the phase diagram of this (or related) ring Hamiltonian. 

The paper is organized as follows. In Scc.|TT]we discuss 
in more detail the wavefunction motivation for consider- 
ing the DBL states. In Sec. IIIII we introduce the lattice 
gauge theory description and study properties of the DBL 
state in the mean field theory that ignores gauge fluc- 
tuations. This provides an initial guide to the singular 
(Bose) surfaces, and is followed in Sec. IIVI with numeri- 
cal characterization of the properties of the actual DBL 
wavefunctions. In Scc.[V] we consider the full gauge the- 
ory description, focusing on the effects of the gauge fluc- 
tuations on the singularities across the Bose surfaces, and 
in particular obtain the long-distance properties such as 
Eq. ((5]) and Eq. JJJ) to order 1/N in our systematic large 
N approach. In Sec. IVII we address the issue of the sta- 
bility of the DBL by considering a putative fixed point 
theory. We conclude in Sec. I Villi with a discussion of 
physical properties and possible future directions. 



II. BOSON WAVEFUNCTIONS AND NODES 

In this section, we expand on the initial motivation for 
considering the DBL states as a way to perform (a kind 
of) flux attachment transformation in a time reversal in- 
variant manner. We also discuss the nodal structure of 
the bosonic wavefunctions, providing some justification 
for the qualifier "D-wave" in the suggested wavefunction 
names in this work. To this end, we consider a "relative 
single particle wavefunction" - more precisely, a cross- 
section of the many-body wavefunction - defined as fol- 
lows. Fixing the positions of all the particles except one, 
we define a function which depends explicitly on the co- 
ordinates of the one "test" particle, and implicitly on the 
coordinates of all the other particles, 

$ r2 ,..., rjv (r) ee V(r,r 2 ,...,r N ) . (11) 

For notational ease we will henceforth drop the implicit 
dependence on the 2(N — 1) spatial coordinates and just 
use the notation $b(r). 

A. Laughlin v = 1/2 revisited 

Consider the Laughlin wavefunction for bosons in a 
half- filled Landau level, Eq. (fT]). In this case, the rela- 
tive single particle wavefunction &(z) is complex and has 
double strength zeroes at the positions of all the other 
-/V — 1 particles. Upon close approach to a specific parti- 
cle, say Zi, the function &(z) ~ (z — Zi) 2 is of a d+id form. 
It is in this sense that the v — 1/2 Laughlin state can 
be viewed as a D-wave fluid: All pairs of particles, upon 
close approach in real space, are in a two-body d + id 
state. Our goal is to construct a time reversal invariant 



analog of the v = 1/2 Laughlin state in which particle 
pairs are similarly in a relative d xy or d x i_ y i state. 

In thinking about Laughlin states, it has been partic- 
ularly instructive to view them in terms of composite 
particles created by "flux attachment." For example, if 
one re-expresses the v = 1/2 Laughlin state as, 

*,=V2 = n> - *i) *c/ , (12) 

i<j 

the "composite fermion" wavefunction is simply a (Van- 
dermonde) determinant, det V, of filled Landau level or- 
bitals, 

# c/ = detV = Hizi-Zj) , (13) 

with Vij — (zi)' J . In a second quantized framework, 
flux attachment can be achieved by introducing a Chern- 
Simons gauge fielcU^. Chern-Simons gauge theory has 
been a useful tool to describe the full Haldane-Halperin 
hierarchy of fractional quantum Hall states, encoding the 
fractional charge and statistics of the quasiparticles as 
well as the structure of the edge excitation a 19 i 20 . But 
in discussing bosons in a time reversal invariant set- 
ting (zero magnetic field) , flux-attachment techniques are 
problematic. For instance, the usual flux smearing mean 
field approximation is likely to inadvertently break time 
reversal invariance. 

In seeking to avoid Chern-Simons theory, it is worth 
noting that the v = 1/2 state is a perfect square - a 
square of the Vandcrmonde determinant, 

*,=i/ 2 = (detV) 2 . (14) 

This suggests that the v = 1/2 Laughlin state can be 
fruitfully described within a slave particle framework^ 
Consider decomposing the boson creation operator as a 
product of two fermions, b^(r) = t4(r)<4(r). As dis- 
cussed in the Introduction, there is a local gauge redun- 
dancy, and a correct treatment requires the presence of 
a gauge field minimally coupled to both fermions. In the 
slave particle mean field approach to the v = 1/2 boson 
problem, one simply drops the gauge field, obtaining a 
problem of two fermion flavors each moving in a mag- 
netic field. If the electrical charge of the boson is divided 
equally, each fermion flavor is effectively in a full Landau 
level. The mean field wavefunction is, 

^/^n^-^xii^-^)- as) 

i<j i<j 

To obtain a wavefunction for the bosons it is necessary 
to project into the physical Hilbert space, and this is 
achieved here by simply setting = = Z{. One 
thereby recovers the Laughlin state as a square of the 
Vandcrmonde determinant, Eq. (|14[) . 
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B. Time Reversal Invariant Wavefunctions 

Within a first quantized framework, the relative sim- 
plicity of a time reversal invariant bosonic superfluid as 
compared to the Laughlin state is the nodelessness of the 
ground state wavcfunction. As Feynman argued many 
years ago, for non-relativistic bosons moving in the con- 
tinuum with an interaction only depending on the parti- 
cle positions, the kinetic energy of any bosonic wavefunc- 
tion which has sign changes could be reduced by mak- 
ing all the signs positive while keeping the magnitude of 
the wavefunction fixed to leave the potential energy un- 
changed. The ground state wavefunction should thus be 
everywhere non-negative. For the non-interacting Bose 
gas the ground state wavefunction is just "J = 1, but in 
the presence of interactions a popular variational wave- 
function is of the Jastrow form^ 

*jastrow(ri,r 2 ,...,r w ) cx e -£i<i , (16) 

with variational freedom in the two-particle pscudopo- 
tential it(r^ — Tj), which is usually taken to approach 
zero as l/|r| p at large separations. 

As in the case with a magnetic field present, we again 
consider the relative single particle wavefunction. For 
a Bosc condensed superfluid in a time reversal invariant 
system $f>(r) can be taken as real and is everywhere non- 
negative. If there are repulsive interactions between the 
bosons in the superfluid, the amplitude of $b(r) will be 
reduced when the test particle is taken nearby another 
particle, which is implemented by the Jastrow pseudopo- 
tcntial in Eq. (jT5J) . But the sign of $b(r) will remain 
positive, so in some sense all of the particle pairs are in 
a relative S-type state. In the special case of a hard core 
interaction, <E>;,(r) — > for r — > r<, so one can then view 
this as an "extended" S-type wavefunction. 

Motivated by the preceding discussion of the v = 1/2 
Laughlin state, in Eq. (jSJ) we introduced a simple time re- 
versal invariant wavefunction for hardcore bosons which 
is the square of a determinant for free fermions fill- 
ing a Fermi sea. By construction this wavefunction is 
non-negative. Moreover, the wavefunction will have ze- 
roes which coincide with the nodes of the filled Fermi 
sea, ^fps(r\, r 2 , . . . , rjv) = dcte lki ' rj . With time re- 
versal invariance, a relative single fermion wavefunction, 
<fr/(r) = ^Fs( r , r 2, ■ ■ ■ ! r iv)j can be taken as real. As a 
result, the nodal structure of 4>/(r) will be qualitatively 
different than for the filled Landau level state, vanishing 
along nodal lines rather than at isolated points. These 
nodal lines will pass through the positions of all of the 
other fermions^ Upon taking the "test" particle across 
a nodal line, the function $/(r) changes sign, vanishing 
linearly upon approaching the nodal line. If one takes 
the position of the test particle close to another particle, 
$/(r) will have a p-wave character - in particular a p x 
form if we define the x-axis as being parallel to the nodal 
line. 

Despite the p-wave character of the free fermion deter- 
minant, the determinant-squared wavefunction will not 



have a d-wave character. Rather, the relative single bo- 
son wavefunction will vanish quadratically upon crossing 
the nodal lines of the fermion determinant. When the 
test particle is taken near another particle r,, the func- 
tion vanishes quadratically, <J>i>(r) ~ (x—Xi) 2 . The "pair" 
wavefunction is thus of an "extended S-type" form, van- 
ishing along the residual nodal line. 

Let us now consider the DBL many-body wavefunction 
Eq. (|5J), which is the product of two different Slater de- 
terminants for fermions that fill elliptical Fermi surfaces 
as shown in Fig. [1] The nodal structure is revealed by 
exploring $b(r). This function will have two sets of nodal 
lines, one from each of the determinants. Due to the el- 
liptical nature of the respective Fermi surfaces, the two 
sets of lines, both of which pass through all of the parti- 
cles, will generally not coincide with one another. Indeed, 
the fermion determinant coming from an elliptical Fermi 
surface will have nodal lines running preferentially per- 
pendicular to the long axis of the ellipse. Focusing on 
the behavior of 4>(,(r) near a target particle r,-, one antic- 
ipates a behavior of the d xy form, $(,(r) ~ (x—Xi)(y—yi). 
Here we have assumed that the two nodal lines actually 
coincide with the x and y axes. In general, for a typi- 
cal configuration of fixed particle coordinates and a given 
target particle, this will not precisely be the case. More 
generically, the two nodal lines will intersect a particle at 
two angles which are not aligned with the axes. But the 
sign of the relative wavefunction when the test particle 

encircles the target particle will still behave as H 1 — , 

the same sign structure as a D xy or D x 2_ y 2 orbital. This 
is illustrated in Fig. O 

Since the DBL many-body wavefunction is not node- 
less, it cannot be the ground state of a continuum Hamil- 
tonian of bosons. If we put the coordinates on the sites 
of a 2d square lattice, the ground state wavefunction is 
only assured to be non-negative if the sign and form of 
the hopping matrix in the lattice tight binding Hamilto- 
nian is such that it satisfies the Marshall sign conditions 
- the requirement that a choice of gauge is possible to 
make all of the off-diagonal matrix elements negative. In 
Sec. I VIII we consider a particular Marshall sign violating 
lattice boson Hamiltonian which might exhibit a ground 
state of the proposed D-wave form. 



C. Precedents of (det)i x (det)2 wavefunctions for 
spin liquids 

Before focusing solely on the DBL states, we want 
to mention that our construction of time reversal in- 
variant bosonic wavefunctions as a product of two dis- 
tinct determinants has nice precedents in the studies 
of spin liquids on the triangular lattice. One can view 
the triangular Heisenberg antiferromagnet as a system of 
hard-core bosons at half-filling in the background field of 
flux 7r through each triangle. Kalmcycr and Laughlin 22 
proposed to view this in the continuum, obtaining a 
boson system at v = 1/2, and their chiral spin liq- 
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uid wavcfunction is precisely the lattice analog of the 
v = 1/2 state Eq. |T]). Alternatively, using a slave 
fermion approach ) 23 ' 24 = d\d\, we divide the boson 
charge equally between the two fermions, so each sees on 
average a flux of 7r/2 per triangle; the mean field where 
the d\ and d 2 see the same static flux of ir/2 per triangle 
gives a filled Landau level for each fermion and repro- 
duces precisely the Laughlin-Kalmeyer chiral spin liquid. 

However, we can be more creative about the fluxes 
seen by the slave fermions while maintaining the average 
of 7r/2 flux per triangle. One example is to take differ- 
ent flux patterns for the two species as follows: For the 
d\ fermions, put flux through all up-pointing triangles 
and 7r flux through all down-pointing triangles, while for 
the d 2 fermions interchange the locations of the and 
7r fluxes. This state is in fact identical to the so-called 
UXBt^t^t^ spin liquid found in Ref . I25I. It is a time re- 
versal invariant gapless algebraic spin liquid (ASL) with 
Dirac nodes in the spinon spectrum; it has very good en- 
ergetics for the nearest neighbor triangular antifcrromag- 
net, starting from the isotropic lattice and all the way to 
the limit of weakly coupled chains (e.g., Ref. [26| found a 
different gauge-equivalent formulation of this state with- 
out realizing its ASL character). 

Another example is obtained by taking yet different 
flux patterns: Select one lattice direction - chain direc- 
tion in the anisotropic lattice case. For the d\ fermions, 
put flux through triangles siding even chains and 7r 
flux through triangles siding odd chains, while for the d 2 
fermions interchange the locations of the and 7r fluxes. 
This is again a time reversal invariant ASL and is listed 
as UICtIt^t 1 in Ref. HI Unlike the UIB state, there is 
no isotropic U1C liquid, but the energetics performance 
of U1C and UIB spin liquids is almost indistinguishable 
for weakly coupled chains and matches that of compet- 
ing magnetically ordered states; in this context, the time 
reversal invariant UIB and U1C liquids are significantly 
better than the chiral Laughlin-Kalmeyer variant. 



Our construction of the DBL states differs in that 
we do not require any special filling for the bosons and 
there is no special Heisenberg spin symmetry. Neither 
of these special conditions of the spin model setting are 
needed for the construction and subsequent gauge theory 
analysis to go through. Given the growing belie f 11 ' 27 ! 28 
that critical spin liquids do exist, we do not see any 
reasons why the situation should be any different for 
our boson liquids at arbitrary incommensurate densities. 
In our construction, such liquids will gencrically have 
some underlying partially filled bands and therefore 
Fermi surfaces of slave fermions. Of course, whether a 
particular state is realized in a given model requires a 
detailed case by case study, and in Sec. IVIII we suggest 
some frustrated boson models that may stabilize the 
DBL phase. Our primary goal in the next Sees. IIIIHVII 
will be to characterize the DBL states without worrying 
where to find them. 

III. MEAN FIELD THEORY FOR THE D-WAVE 
BOSE LIQUID 

A. Gauge Theory formulation 

We next consider the challenge of constructing a field 
theory which can access such a .D^-Bose liquid (DBL) 
state. Due to our inability to implement flux attach- 
ment in a tractable time reversal invariant manner, we 
follow instead the slave particle approach. As above, we 
decompose the hard core lattice boson as a product of 
two fermions, = d\d\. Consider then a lattice U(l) 
gauge theory on the square lattice in terms of these slave 
particles: 

Hu(i) — H t + H a , (17) 
with the fermion hopping Hamiltonian of the form, 



J 



H t = - 



[*||e w * (r) (4(r)di(r + x) + t ± e ia » (r) 4(r)di(r + y) + h.c. 

r 

[txe- ia *( r) dl(r)d 2 (r + x) + t l{ e~ ia ^ d\{r)d 2 {v + y) + h.c. 



(18) 
(19) 



The d\ fermion hopping amplitudes in the x and y direc- 
tions are t\\ and t± respectively, while the two amplitudes 
are interchanged for the d 2 fermions. In the following, 
we take ty > t± for concreteness. Note also that the two 
fermion species carry opposite gauge charges. The gauge 
field Hamiltonian is simply 



ffo^EE e » - ^'Ecos[(V x a), 



(20) 



r fi—x,y 



where the lattice "magnetic" field is 

(V x a) T = a x (r) + a y (r + x) - a x (r + y) - o„(r) . (21) 

The integer- valued "electric" field e M (r) is canonically 
conjugate to the compact gauge field a Al (r) on the same 
link. The above Hamiltonian is supplemented by a gauge 
constraint on the physical states, which must satisfy 
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Gauss' law, 



(V-e) r = 4(r)di(r)-4(r)d2(r) , 



(22) 



with 



(V ■ e) r = e x (r) - e x (r - x) + e. y (r) - e y (r - y) . (23) 

In the limit h K,t\\,t±, the electric field vanishes, 

and Gauss' law reduces to d\d\ = d\di which projects 
back into the physical boson Hilbert space. In this strong 
coupling limit, it is possible to perturbativcly eliminate 
the gauge field to obtain a Hamiltonian for hard core 
bosons hopping on the square lattice with additional ring 
exchange terms. We pursue this in Section [VIII where we 
compare the energetics of the DBL wavefunction with 
other states. Here we instead focus on the weak coupling 
limit, K ^> h, which suppresses the magnetic flux. The 
usual slave particle mean field treatment corresponds to 
simply setting (V x a) r equal to a constant. The sim- 
plest mean field state, and the one which should corre- 
spond to the wavefunctions in the previous section, is 
with zero flux through all plaquettes. The corresponding 
mean field Hamiltonian describes non-interacting slave 
fcrmions. Each species has anisotropic near-neighbor 
hopping amplitudes, but the two are related under the 
90 degree rotation, thus producing a boson liquid that 
respects the symmetries of the square lattice. 



B. Mean field results for the DBL 

In order to focus on the effects of the underlying Fermi 
surfaces without the complications of lattice physics such 
as Brillouin zone folding, we first consider fermions in 
the 2d continuum with anisotropic effective masses. The 
resulting elliptical Fermi surfaces are, 



di : {wk x ) 2 + {ky/w) 2 
d 2 ■■ {k x /w) 2 + (wk y ) 2 



k 2 

k 2 



(24) 
(25) 



where the parameter w > 1 characterizes the degree of 
eccentricity (the conventional eccentricity of the ellipses 
is given by e = y^l — 1/w 2 ). In the D-wave Bose liquid, 
if signifies the mismatch between the two Fermi surfaces, 
and we will refer to this measure as "D-eccentricity" . We 
study equal time correlation functions since these can be 
compared directly with the properties of the wavefunc- 
tions, which is done in Sec. lIVl we also consider temporal 
dependencies within the mean field as a measure of spec- 
tral properties. It is useful to have in mind that much 
of the following analysis of long-distance properties needs 
only the knowledge of relevant Fermi surface patches and 
not of the full surfaces. 

Consider first the one-particle off-diagonal density ma- 
trix (or Green's function) for the boson, Gf,(r, r) = 
(6^(r, t)£>(0, t)), with Gb(0, 0) = p - the average boson 



density. It will also be of interest to consider the momen- 
tum occupation probability, 



(n k ) = (b{b k ) = / dr G b {v)e 



(26) 



Within the mean field theory, the natural approximation 
for the order parameter correlation is 



G¥ F ( 



G% f {t,t)G% f (t,t)/p, (27) 



where G^ IF are the mean field (bare) fermion Green's 
functions. This approximation satisfies G&(0) = 
Gd a (0) = p. Here and below, the imaginary time r is 
understood to be zero if not explicitly present. 

The fermion Green's functions are readily calculated. 
Thus, at long distances r 3> hp 1 , the main contribution 
to G^ IF (r) comes from the Fermi surface patches where 
the group velocity is parallel or antiparallel to the obser- 
vation direction f = r/|r|. With inversion symmetry, we 
can denote the corresponding patch locations as ±k Fa 
and the Fermi surface curvature as c a , and obtain 



g¥ f (v) 



1 



cos(k/ 



3tt/4) 



2 1 /2 7r 3/2 



c„ /2 |if/ 2 



(28) 



It is important to remember that k Fl , k F2 , and c\ , c 2 de- 
pend implicitly on the direction f and are different in 
general when the d\ and d 2 Fermi surfaces do not coin- 
cide. 

For the equal-time boson Green's function we thus get 



cos[(k Fl - k F2 ) ■ r] 

cos[(k Fl + kp 2 ) • r — 37r/2] 
c\ ,2 c\ ,2 \r\3 



(29) 



which decays algebraically while oscillating with r- 
direction-dependent wave vectors k Fl + kp 2 and k Fl — 
kp 2 . Such wave vectors, which are constructed by consid- 
ering patches on the two Fermi surfaces that are parallel 
to each other, will span some new loci in the momentum 
space as illustrated in Fig. [3] The above large-distance 
behavior corresponds to singularities ~ |i5fc| 3//2 across 
these lines. In the zero eccentricity limit (i.e., when 
w = 1), the two Fermi surfaces coincide, giving 



r<MF 



(r) 



1 + cos[2fc F r - 3tt/2] 



(S-typc) . (30) 



Compared to the general case with w > 1, the k Fl — k F , 
locus shrinks here to zero momentum, while the singu- 
larity in rik is hardened to |fc|. 

As we show in Sec. |V] the power law decay of the mean 
field boson Green's function survives in the presence of 
gauge fluctuations (but with modified exponents). The 
algebraic decay is a consequence (and a measurable in- 
dication) of the "criticality" of the DBL and its gapless 
excitations. Of course, by the very construction we ex- 
pect many gapless excitations because of the underlying 
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Fermi surfaces. Some measure of the low-energy spec- 
trum is contained in the time dependence of the local 
boson Green's function, 



G 



(31) 



where vq is the density of states at the Fermi energy 
for each species. The corresponding local boson spectral 
function is 



A^ F (v = 0,E) 



aMF 
A b 



(k,E) 



(32) 



In the discussion of the DBL phase on the lattice in 
Sec. 1111 CI we will encounter a situation where the boson 
Green's function decays exponentially in space because 
of the topology of the Fermi surfaces, while the above 
spectral signature of the low-energy excitations depends 
only on the density of states which is a property of the 
entire Fermi surface and is insensitive to the topology 
otherwise. 

The most natural instability of the DBL is towards 
a superfluid state. As we discuss in Section IVI1 in the 
limit of vanishing D-eccentricity (w = 1), the resulting 
S-type Bosc Liquid phase obtained in mean field theory 
is, in the presence of gauge fluctuations, most probably 
unstable to superfluidity. Moreover, our Sec. IIVI anal- 
ysis of the determinant squared wavefunction appropri- 
ate to the S-Bose Liquid is consistent with off-diagonal 
long-range order. These results strongly suggest that the 
mean field S-type Bose liquid phase probably can not 
exist as a stable quantum phase. However, with non- 
vanishing D-eccentricity, both the gauge theory analy- 
sis in Sec. I VII and the properties of the corresponding 
(det)a;(det)j, wavefunction which we explore in Sec. IIVI 
suggest that the DBL is a stable critical quantum phase. 

It is also instructive to examine several other correla- 
tion functions within the present mean field treatment. 
Specifically, consider the boson density-density correla- 
tion function, 



D b (r) = (:6p(p)6K*)--) 

= N{N-l)(S d {r 1 )6 d {r 2 



(33) 
(34) 



where N is the total number of particles. As defined, 
D b (v) approaches zero for large separations, while nega- 
tive values at small distances signify a correlation hole. 
The density structure factor can be calculated as 



D b (k) = / dr D b (r)e 



-ikr 



v 



p 



(35) 



where 8pu = exp[— ik-Tj] — NS^^q and V is the system 
volume. 

Microscopically, since W = d\d\, for each boson added 
to the system both a d\ and a d 2 fermions are added. 
As such, the boson and fermion densities are equal: 
Pb = Pdi = Pd-2- However, in the mean field treatment, 
the two fermion flavors propagate independently with 




k x /k F 

FIG. 3: Singular lines in the momentum space for the boson 
Green's function (left) and the density correlation (right). In 
each panel, the thin lines show the elliptical di and di Fermi 
surfaces for the parameter w 4 = 3, cf. Eq. (|24|l . In the left 
panel, the k.F 1 ikp, loci are constructed by considering a pair 
of di and di fermions with parallel (+ sign) or antiparallel ( — 
sign) group velocities. The axes are scaled by an "average" 
kF that would obtain if the Fermi surfaces were circular for 
the same particle density. 



different Fermi surfaces, so the corresponding fermion 
density-density correlation functions, which we denote as 
D^ IF , do not coincide. This ambiguity in the density cor- 
relation is an intrinsic deficiency of the mean field theory. 
As a crude measure, we approximate the boson density 
correlation function as an average of that of the d\ and 
(f 2 fermions, 



(36) 



The individual density correlation functions in the mean 
field theory are simply, 



n MF 



(r) 



l + cos[2k Fa -r-37r/2] 



(37) 

where again kp a is the place on the d a Fermi surface 
where the normal is parallel to the observation direction 
f . We recognize the 2kp oscillation with power law en- 
velope, which in the momentum space translates to a 
singularity in the structure factor ~ \2kp — fc| 3 / 2 across 
the 2kp line (such 2kp surfaces are illustrated in Fig. [3]), 
while there is also a \k\ singularity at zero momentum. 
As obtained from Eq. (|36[) . the mean field boson density 
correlator D^ IF has singularities at both 2kp ± and 2kp 2 , 
and "knows" about the presence of both Fermi surfaces. 
Remarkably, this mean field approximation is consistent 
with our analysis of the (det) 2 -(det). y boson wavefunction 
in Sec. IIVI which also reveals singular behavior at both 
2kF 1 and 2kp 2 . 

Let us now consider the two-boson correlator. 



GafcCn^ri.r'a) = (&t(r 1 )fet(r 2 )6(r' 1 )&(r' 2 )> 



(38) 



which injects a pair of bosons at ri , r 2 and removes a pair 
at , r 2 . Consider first the limit where the separation be- 
tween the two injected bosons, r = ri — r 2 , and the two 
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removed bosons, r' = r[ — r' 2 , are both small compared 
to the mean inter-particle spacing, |r|, |r'| <C kp . More- 
over, take the separation R = ri j^ r2 — between the 
injected and removed pair to be much larger than the 
inter-particle spacing, |R| 3> kp 1 . In this limit the mean 
field 2-boson correlator can be expressed as, 

c mf k F ^ y (r;R)^„(r^;R) 

^ ~ [l+W sin 2( 2 ^)]5/2 (k F R)6 >W 

with "pair wavefunction" , 

* KJ/ (r; R) = (r ■ R) 2 + W sin(20)a;y . (40) 

Here for concreteness we specialized to the elliptical 
Fermi surfaces, Eq. and W > characterizes the 

D-eccentricity, 

W= wA + ^-l; (41) 

R = R/i? = cos(0)x + sin(</>)y, is the unit vector 
pointing from one pair to the other. With vanishing 
D-eccentricity, W = 0, the pair wavefunction is of an 
"extended-s" form, for example ^> xy ~ x 2 for = 0. This 
corresponds to a quadratic nodal line in the pair wave- 
function. With W = the directions of the nodal lines 
for each pair are aligned perpendicular to the vector con- 
necting one pair to the other. On the other hand, in the 
limit of very large D-eccentricity, W 3> 1, the pair wave- 
function takes the d xy form, ^ xy ~ xy, corresponding to 
two nodal lines along the x and y axes. Upon rotating 
the internal coordinate, the sign of the pair wavefunction 
takes the usual d xy form, precisely as for the Cooper pair 
wavefunction in a d xy superconductor. Thus, the DBL 
has quasi-long-range order in the d xy boson pair chan- 
nel, although the power law decay exponent within mean 
field theory is very large. 

It is also interesting to examine the time decay of a 
"local" two-boson Green's function, 

G 2b (v, r'; t-t')^ (b\r, rp (0, r)b(v', r>(0, r')) , 

(42) 

with |r|, |r'| < kp. A pair of bosons is injected at posi- 
tions r and and is removed at later times at positions 
r' and 0. Within mean field theory, 

G^(r, r'; r) ~ l[(r • r') 2 + 2W(xy)(x'y')] . (43) 

For large D-eccentricity, W 3> 1, this factorizes into a 
product of two "pair wavefunctions" of the d xy form. 
This correlator can be used to extract a local pair tun- 
neling density of states, 

Pxy {E) = J (Vt y (T)V xy (0))e- ElTl dT , (44) 
where ~D xy injects a d xy pair centered at the origin, 

t>i y = [ xyeM-\r\/tiW(r)b\0) , (45) 



with pair "size" £. Within mean field theory one obtains 
a power law tunneling density of states p xy F (E) ~ WE 3 , 
with an amplitude that grows with the D-eccentricity pa- 
rameter W. The tunneling density of states for an s-wave 
pair also vanishes as E 3 but the amplitude is independent 
of W. This pair tunneling density of states is perhaps the 
best diagnostic for measuring the degree of local d-wave 
two-boson correlations in the DBL. 

Finally, we consider a box correlator for bosons which 
is defined in the Introduction, Eq. © . Within mean field 
this factorizes as, B^ IF (x) = Bd 1 {x)Bd 2 {x), and from 
Wick's theorem and 90-degree rotational invariance, 

B dl (x) = G\ (0, x) - G\ (x, 0) = -B d , (x) . (46) 

Thus, 

B£ IF {x) = -[Gl(0,x) - G 2 di (x,0)} 2 < • (47) 

Notice that in the case with zero eccentricity the mean 
field box correlator vanishes, but since the (det) 2 wave- 
function is non-negative, one expects Bb to be positive 
upon inclusion of gauge fluctuations. With non-zero D- 
eccentricity, the mean field box correlator is negative 
for all spatial separations, which reflects the underlying 
d xy nodal structure, and decays as x~ e . It is possible 
that the inclusion of gauge fluctuations will again mod- 
ify the mean field result in the case with closed Fermi 
surfaces because of the pairing tendencies of the d\ and 
di fermions. In the case with open Fermi surfaces to be 
discussed next, we conjecture that the exact box corre- 
lator will be negative at large spatial separations while 
decaying to zero as the box size is taken to infinity. This 
conjecture appears to be consistent with the box corre- 
lator extracted from the (det) a; (det) !/ wavefunction with 
open Fermi surfaces in Sec. IIVB1 

C. Case with open Fermi surfaces 

The preceding analysis holds for arbitrary Fermi sur- 
faces, in particular, for the lattice bands obtained from 
Eq. (fT5| . We only need to remember that the fermion 
Green's function Eq. (|28|) is determined by the Fermi 
surface patches with the group velocities that are paral- 
lel or antiparallel to f . Once Gd t and Gd 2 are known, the 
other correlation functions follow. 

An interesting situation occurs on the lattice when the 
ratio in /t± is such that the d\ and d^ Fermi surfaces are 
open, which is illustrated in the right panel of Fig. [TJ In 
this case, for an observation direction close to, say, the 
y-axis, there are no d\ Fermi surface patches with nor- 
mals in this direction, so the d\ Green's function has an 
exponential decay in this direction instead of the power 
law Eq. Since the real-space boson Green's func- 

tion is the product of the two fermion Green's functions, 
Eq. (|2"T|) , we conclude that it decays exponentially in the 
directions near the x- and y-axes. There may still still be 
directions near x ± y in which both Gd 1 and Gd 2 show 
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power law behavior, and so will the boson Green's func- 
tion. 

Eventually, for large enough t\\/t±, the two Fermi sur- 
faces will have no parallel patches, and the boson Green's 
function will decay exponentially in all directions. We 
will call this phase DLBL for D-wavc correlated Local 
Boson Liquid. Note, however, that the system is still 
gapless and critical, as can be measured, e.g., from the 
local spectral function Eq. (f3"2")h Also, the boson density 
correlations still have the power law envelope, Eq. (f37|) . 
if one fermion field can propagate in the observation di- 
rection. Furthermore, the boson box correlator, Eq. (|4"T)l . 
exhibits power law behavior even though the single and 
pair boson Green's functions arc exponentially decaying. 
From the boson field perspective, the system is local and 
bosons have hard time to propagate; nevertheless, the 
system has power law correlations in other properties 
and, in particular, charges can propagate. 

Based on an analysis of gauge fluctuations in Sec. [Vj 
we conjecture that the d\ and d,2 systems effectively de- 
couple at low energies in the DLBL, and at large dis- 
tances the exact box correlator is negative and decays to 
zero as a non-oscillatory power law with an integer expo- 
nent which is independent of non-universal lattice scale 
physics, Bb(x) ~ — x~ s . This is the same as in the mean 
field theory except with a larger power. 

It is also worth mentioning the limiting case t± = 
that gives completely flat Fermi surfaces, which we will 
call "extremal DLBL." The d\ fcrmions can move only 
along the x-axis, 

„ , „ sinffcpx) , , . 

G dl (x,y) = S y0 — , (extremal DLBL) , (48) 

TTX 

while the c?2 fermions can move only along the y-axis. As 
a result, the boson field cannot propagate at all, not even 
one lattice spacing. However, this special system still has 
power law density-density correlations, e.g., 

D% F (x, y) ~ -fl„ |0 gg-ffjgg) , (extremal DLBL) , 

(49) 

as well as power law box correlation, 

B A b IF (x) ~ - Sin4( ^ ) , (extremal DLBL) . (50) 

a; 4 

From the numerical study of the extremal DLBL wave- 
function, Sec. lIVB~2l we conjecture that these mean field 
power laws also hold upon the Gutzwillcr projection. In 
the gauge theory context, we would say that the fermions 
remain unaffected by the gauge field fluctuations. This 
extremal DLBL wavefunction is of interest because of its 
similarity to the so-called Excitonic Bosc Liquid ground 
state of a pure ring exchange model - we discuss this in 
Sec. EE 



IV. PROPERTIES OF THE WAVEFUNCTIONS 

In this section we study the \& = (det) x (det) 2/ wave- 
functions directly and measure numerically their proper- 
ties such as the boson Green's function and the density 
correlation defined in Sec. IIII Bl A detailed comparison 
is made with the mean field, which provides an initial 
guide. This is followed by interpretations of the observed 
deviations from the mean field using simple "Amperean 
interaction" rules of thumb; the actual calculations be- 
hind these rules within the gauge fluctuations theory are 
given in Sec. fVl 

We first consider wavefunctions in the continuum so as 
to avoid lattice effects and focus on the consequences of 
the underlying Fermi surfaces. It appears that the S-type 
wavefunction, = (det) 2 , has off-diagonal long-range or- 
der, while a generic DBL wavefunction with non-zero D- 
cccentricity has only power law boson correlations. We 
then consider wavefunctions on the lattice where we can 
access the DLBL phase with open Fermi surfaces de- 
scribed in Sec. IIII C) this features boson Green's func- 
tion that decays exponentially in real space. Finally, we 
discuss in some detail the extremal wavefunction that is 
obtained when the Fermi surfaces are completely flat. 

The calculations with the wavefunctions are performed 
numerically using standard dctcrminantal VMC (varia- 
tional Monte Carlo) techniques^ The system used in all 
calculations is a square box with periodic boundary con- 
ditions. In each case, the particle number is chosen so as 
to fill complete momentum shells under the Fermi sur- 
faces. To facilitate the comparison, the presented mean 
field is calculated for the same finite systems. 

A. Wavefunctions in the continuum 

Before proceeding with the numerics, the DBL struc- 
ture factors can be found analytically in some ranges in 
the momentum space: Thus, by expanding the boson 
wavefunction in terms of the orbitals that form each de- 
terminant, it is easy to see that vanishes outside the 
kpi + surface of Fig. [31 while D(k) vanishes outside 
the 2kpi + 2ki?2 surface. 

1. S-type state 

We begin by considering the limit with zero D- 
cccentricity, in which case the boson wavefunction is pos- 
itive everywhere except for the nodes. Figure [¥] shows 
boson Green's function measured for two systems with 
N = 161 and 325 particles. It appears that at large sep- 
arations Gb(r) approaches a finite positive value which 
is roughly the same for both system sizes. To be more 
quantitative, the k = mode contains n^o/N = 0.10 
and 0.083 fractions of bosons in the two systems. The 
present data extracted from the continuum wavefunction 
cannot rule out the possibility that the Green's function 
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S-type wavefunction, *P = (det) 2 
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x / a 



FIG. 4: Boson correlation in the (det) 2 wavefunction for two 
systems with N = 161 and 325 particles. The distance is 
measured in units of a = / o _1//2 , and the plots are cut when 
the separation exceeds half of the box size. Since Gb(r) main- 
tains a sizable limiting value that depends only weakly on 
N, we suggest that there is off-diagonal long-range order in 
this wavefunction. The mean field result is also plotted, and 
the vertical scale is chosen to accommodate its fast decay, 
G b ~ l/r . 



vanishes in the large distance limit. But at the very least, 
the result of the projection is clearly dramatic in this 
case with matched Fermi surfaces, since the fall-off of 
the boson correlation is very slow (if any). We also re- 
mark that we unambiguously find off-diagonal long-range 
order when this wavefunction is studied on a lattice at 
fixed boson density 

The extremely strong enhancement of the boson cor- 
relation over the mean field prediction seen in Fig. [4] can 
be understood qualitatively as a result of the pairing of 
the d\ and di fermions mediated by the gauge field. In- 
deed, as we discuss in Sec. [V] the constituents of such 
a (zero momentum) "Cooper pair" with oppositely di- 
rected group velocities and opposite gauge charges pro- 
duce parallel gauge currents and therefore experience 
Ampcrean attraction mediated by the gauge field. This 
"Amperean attraction" rule of thumb for the enhance- 
ments in correlations appears to be taken to the extreme 
in the (det) 2 wavefunction, where we can crudely picture 
the fermions paired back into b and condensed, giving 
rise to long-range order in Gb- We also point out that 
this wavefunction has rather unusual density-density cor- 
relations which arc singular around the 2kp circle. Since 
the presence of the off-diagonal order makes this state 
less interesting to us, we do not consider such details any 
further. The observation of the order suggests that the 
S-type state is unstable; the ultimate phase in this case 
is likely a conventional superfluid, and a good superfluid 
wavefunction needs to be constructed differently, e.g., us- 
ing Jastrow approach that builds in proper density cor- 
relations. 



DBL wavefunction, *P = (det) x (detL 

^ D(k) 
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FIG. 5: Top: Mode occupation Uk and density structure 
factor D(k) measured in the continuum DBL wavefunction 
with w 4 — 3 (note that our -D(k) defined in Eq. I35I vanishes 
for large k). The boson number is N = 325; the underlying 
Fermi surfaces and all singular lines are shown to scale in 
Fig- Bottom: Cut through the data in the (1,0) direction 
in the k-space and comparison with the mean field results. 
In the left panel, the arrows point the appropriate kpi ± fcj?2 
momenta (see Fig. [3]), while in the right panel the arrows 
point the 2/cfi and 2&F2 momenta. 



2. DBL state 

We now consider in detail a representative case with 
non-zero D-eccentricity - specifically, with w 4 = 3 in 
Eq. (|24|) . The corresponding d\ and ^2 Fermi surfaces are 
shown to scale in Fig. [3] together with the singular lines 
for the order parameter and density correlations identi- 
fied in the mean field, Sec. IIIIl 

Figure [5] gives an overall view of and D(k) in the 
two-dimensional momentum space, and also shows a one- 
dimensional cut in the (1,0) k-space direction together 
with the mean field predictions. 

Consider first the mode occupation function n^. Upon 
projection, this develops a squarish top with somewhat 
sharper edges as compared with the mean field. The 
latter is not shown in the full k-spacc but looks more 
smooth, while a (1,0) cut can be seen in the bottom panel 
of Fig. [5] Coming from large momenta, significant devi- 
ations from the mean field set in near the kpi — ki?2 line 
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FIG. 6: Boson Green's function Gb(r) and density correla- 
tion Db{r) measured along the r = (x,0) direction for the 
DBL wavefunction with N = 325 bosons and D-eccentricity 
characterized by w 4 = 3 (this is the same system as in Fig. [5| ■ 
The distance is measured in units of a — p~ x ^ 2 , and the plots 
are cut roughly where the signal in Gb(r) approaches the noise 
level. Note the logarithmic scales used, so it is the absolute 
values that are plotted, while the positive/negative signs of 
the data are indicated with filled/open symbols respectively. 
In the case of the mean field data, the signs are indicated us- 



ing solid/broken lines; D is always negative, see Eq. ([37 



of Fig. [3] (the corresponding location in the ID cut is in- 
dicated with an arrow). This surface is singular in the 
mean field, but the singularities are weak and are almost 
not visible, while they become more pronounced upon 
the projection. On the other hand, we observe no such 
enhancement near the kpi +i^F2 line. 

This difference between the kpi — k^2 and kpi +kp2 is 
also visible when examining the boson Green's function 
in real space as shown in the left panels of figures [5] and [7] 
The mean field Gf F (r), Eq. (2SJ), has both k F i-k F2 and 
kpi + kp2 components equally present, while after the 
Gutzwiller projection one finds that the kpi — k^?2 com- 
ponent dominates. Indeed, for the boson Green's func- 
tion measured along the x-axis, Fig. [51 the relevant Fermi 
surface patches have normals in the (x,0) direction and 
are easily located (the relevant momentum space cut of 
nk is also shown in the left panel of Fig. [5]) . By comparing 
with the mean field Eq. I|29p . one finds that after the pro- 
jection the amplitude of the oscillation with the smaller 
wave vector \1cfi — wins over that with the larger 
wave vector fcpi + fcp2- Slightly more care is needed when 
interpreting Gt,(r) in the x — y diagonal direction, Fig. [7] 
In this case, the appropriate locations on the kpi — kp2 
surface are at the "cusp" points with k x = —k y in Fig. [31 
since this is where the normal to the singular surface is 
parallel to the observation direction f . This component 
does not oscillate when moving along the x — y diagonal 
in real space, and when it is enhanced, the oscillations 
due to the other kpi + kj?2 component become less visi- 
ble, which is what what we see in Fig. [7J 
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FIG. 7: The same as in Fig. [6] but in the x = y diagonal 
direction in real space. 



Unfortunately, from this data we cannot attest whether 
the enhanced singularities are characterized by new ex- 
ponents or whether we see only an amplitude effect. The 
locations where the enhancements occur are consistent 
with the Ampcrcan rules of thumb, Sec. [Vl Indeed, the 
d\ and e?2 constituents of the boson operator = c^cd, 
move in the opposite directions when contributing to the 
correlation at kpi — kj?2 and therefore experience Am- 
perean attraction and enhancement (remembering that 
d\ and di carry opposite gauge charges), while they move 
in the same direction when contributing at k^i+k^ 2 and 
therefore experience Amperean suppression. These pic- 
tures are made more precise in Sec. [VI where we find 
that within the gauge fluctuations theory, the enhanced 
correlations are characterized by new exponents that de- 
pend on the observation direction because of the varying 
degree of the Fermi surface curvature matching in the 
d\-&2 pairing channel. The strongest such enhancement 
is expected along the diagonals x = ±y, which is roughly 
what we find in the projected wavefunctions (we repeat 
again that we cannot make any statement about the ex- 
ponents from the data other than that we see increased 
numerical correlations compared with the mean field). 

Consider now the density correlation D^{r) and the 
corresponding structure factor _D(,(k). In Fig. [5) we can 
see the singular 2kpi and 2k^?2 lines, cf. Fig. [3j the en- 
hancements are peaked where the two curves cross. The 
singular 2fci? points are also visible in the (1,0) k-space 
cut, bottom right panel of Fig. On the other hand, the 
\k\ singularity near zero momentum is not enhanced. Ex- 
amination of the density correlation in real space, right 
panels of figures [S] and [JJ shows an overall increase and 
dominance of the oscillatory components over the zero- 
momentum component as compared with the mean field 
Eq. (|37|) . Again, we cannot tell whether there is a change 
in the exponents or just an amplitude effect. The 2fcp en- 
hancements of the density correlation agree qualitatively 
with the gauge theory expectations. Indeed, consider 
the density operator d\d\. The particle and hole con- 
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stituents (which are oppositely gauge charged) have an- 
tiparallel group velocities when contributing to the 2kp 
component and therefore experience Amperean attrac- 
tion, while they move in the same direction when con- 
tributing at zero momentum and therefore repel each 
other. Again, these rules of thumb are made more precise 
in Sec. [Vj where such enhancements in the particle-hole 
channel are characterized by new direction-dependent ex- 
ponents (the density correlations are found to decay more 
slowly when less curved patches are involved). Note that 
the projection imposes pf, = p<n = Pd 2 i an d we expect 
Db to acquire both the 2fcpi and 2kp2 signatures (in the 
gauge theory, p^ and pd 2 imprint on each other via non- 
singular high-energy connections that are not manifest in 
our low-energy effective description of Sec. Iv|) . 

To summarize, the DBL wavefunction clearly knows 
about the underlying Fermi surfaces; it also contains 
germs of the gauge fluctuations theory, since the enhance- 
ments/suppressions of the various "2fcp" lines appear to 
agree with the Amperean rules. However, there are no 
reasons to believe that the wavefunction and the gauge 
theory will have the same long-distance properties. In 
particular, our measurements can not tell whether the 
long-distance power laws are changed upon the projec- 
tion compared with the mean field. Still, it is gratifying 
to see that the Amperean rules work for the DBL wave- 
function. 



B. Wavefunctions on the lattice 

We now describe our results for various states on the 
lattice. As we have already mentioned, S-type wavefunc- 
tion, = (det) 2 , with finite boson density per site, have 
off-diagonal long-range order, which we confirm unam- 
biguously using finite-size scaling. On the other hand, 
wavefunctions with non-zero D-eccentricity do not show 
such order. While we have not performed detailed stud- 
ies, we expect that the case with closed Fermi surfaces is 
similar to the already described continuum DBL wave- 
function, with a slight complication that one needs to do 
proper Brillouin zone folding when considering singular 
surfaces in the momentum space. 



1. Open Fermi surfaces - DLBL state 

Of particular interest is the case with open Fermi sur- 
faces, which can only be realized on the lattice. Specifi- 
cally, we studied a 24 x 24 system with N = 216 bosons 
(p = 0.375 per site) and the fermion hopping param- 
eters t|| = 1.35, t± = 0.65. In this case, the d\ and 
d,2 Fermi surfaces have no parallel patches; we therefore 
expect that the boson Green's function decays exponen- 
tially. When we measure G&(r), we find that it drops 
below our noise level already at 3 lattice spacings, so the 
corresponding plots are not particularly informative and 
are not shown here. Where we can measure reliably, the 



values after the projection are of the same order as the 
mean field values in the same system. Since the latter 
decay exponentially at large distances, we conjecture the 
same behavior in the DLBL wavefunction. 

On the other hand, the boson density correlation along 
the x and y axes shows oscillations with a clear power 
law envelope, while the correlations are much smaller in 
the diagonal x = ±y directions; this behavior is again 
qualitatively consistent with the mean field, since neither 
d\ nor di Fermi surfaces have normals in the diagonal 
directions, while one or the other has normals along the 
x or y axis. 



2. Flat Fermi surfaces - extremal DLBL state 

Finally, let us discuss the extremal case when the Fermi 
surfaces are completely flat. The d\ fermions can move 
only along the x-axis and di fermions only along the y- 
axis. Using fermion orbitals that are localized on indi- 
vidual rows for d\ (or columns for d%), one can see that 
the boson wavefunction is nonzero only when the number 
of particles on each row (or on each column) is the same. 
Thus, the bosons can not propagate and their inter-site 
correlation is identically zero. As we discuss in Sec. IVII1 
pure ring Hamiltonian conserves boson number on each 
row/column, and the extremal DLBL wavefunction may 
be useful in this context. 

We can still use the box correlator, Eq. ©, to charac- 
terize the state and see some "gaplessness" in the system. 
The measurement is shown in Fig. [51 where we also plot 
renormalized mean field result. The latter is obtained by 
dividing Eq. (|4"T)l by p 2 (l — p) 2 , and a crude justification 
for such procedure^ is as follows: Each d\ or di mean 
field box calculation contains implicitly a weight of or- 
der p 2 {\ — - p) 2 , since for the ring operator to be nonzero, 
two specific sites need to be occupied and two need to be 
empty. However, after the projection it is enough to re- 
quire that only the d\ configuration is "correct" since the 
c?2 fermions are tied to d\. From Fig. [5J we see that the 
box correlator is negative and tracks closely the renor- 
malized mean field values. 

We can also characterize the extremal DLBL state by 
measuring the density correlations and comparing with 
the mean field prediction, Eq. (|4*§)) . The results are in 
Fig. [51 The left panel shows Db(x,0) as a function of 
the distance along the x-axis; it has stronger oscillatory 
component than the mean field and swings back and forth 
across the zero line while the mean field only touches it, 
but the overall magnitudes are comparable and decay as 
1/x 2 . We also find that the density correlation in the x = 
±y diagonal directions decays exponentially (not shown); 
the mean field predicts zero correlation unless strictly 
along the axes, and we expect that after the projection 
this corresponds to exponential decay. 

In the right panel of Fig. [21 we show the density struc- 
ture factor -D(k). Several features are clearly visible: ris- 
ing "towers" draw attention to the lines (±2/cp, k y ) and 
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Extremal DLBL wavefunction; box correlator 




FIG. 8: Box correlator Bt(x) in the extremal DLBL wave- 
function with flat Fermi surfaces on a 24 x 24 lattice with 
TV = 216 bosons (the lines are guide to the eye). Note that 
the box correlator is negative in the mean field and is found 
to be negative in the DLBL state; we plot — Bb(x) to be able 
to use the logarithmic scale. The mean field Eq. (|47|l is renor- 
malized as explained in the text using a numerical factor that 
removes double weighting of configurations. The measured 
values below 10~ 5 are at the noise threshold. 
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FIG. 9: Boson density correlation for the same extremal 
DLBL system as in Fig. [8] Left panel: Real space depen- 
dence along the rr-axis, Db(x, 0), together with the mean field 
expectation. Right panel: Density structure factor in the mo- 
mentum space. 



(k Xi ±2kp); one can also see a "cross" formed by the lines 
(0, k y ) and {k x , 0) that run along the axes. As far as we 
can say, the character of the singularities across these 
lines remains the same as in the mean field; there is an 
amplitude enhancement of the 2kp, but no qualitative 
difference otherwise. In particular, near the line (0,k y ), 
we observe Db, s ing(k x — > 0, k y ) = A\k x \ with A which is 
independent of k y as long as \k y \ ^> \k x \. 

We thus conjecture that the extremal DLBL wavefunc- 
tion is adequately described by the free fermion mean 
field. Some of our findings, e.g., the cross singularity that 



has a long wavelength character, can be understood semi- 
analytically, since the absolute value of the wavefunc- 
tion has a Jastrow form with a peculiar pseudo-potential, 
u(x,y) ~ — ]n\sm(nx/L)\6 y fl — ki\ain(ny/L)\6 Xi o, it 
seems plausible that one can calculate other properties 
of this wavefunction analytically. It is interesting to note 
that the above suggests a free fermion description of this 
boson state, perhaps with some constraints that become 
irrelevant in an infinite system. We also note that the ex- 
tremal DLBL appears to be a relative of the so-called Ex- 
citonic Bose Liquid (EBL) phase predicted in a pure ring 
model on the square lattice^ that we discuss in Sec. I VIII 
the above therefore suggests that there may be some such 
description of the EBL in terms of fermions carrying frac- 
tional 1/2 charge. 



V. GAUGE FLUCTUATIONS 

We now study the gauge theory using analytic tech- 
niques, focusing on the large K limit where it is rea- 
sonable to ignore the presence of magnetic monopoles 
in space-time and to expand the cosine term in Eq. (f^Uj) 
treating the gauge field as non-compact. Moreover, in the 
K — > oo limit, the fermions become free and one can put 
them into a Fermi sea state. A Fermi surface of fermions 
minimally coupled to a non-compact U(l) gauge field has 
a considerable history£2I^^2^2I^^Ml The 
2 + ID such system has been studied most notably as 
a theory for the spin sector in the uniform RVB phase in 
the slave boson approach to the high-T c superconductors. 
It has been argued that such fermion systems have a sta- 
ble phase that in some crude aspects is similar to a Fermi 
liquid - for example, it has a finite long- wavelength com- 
pressibility and spin susceptibility. However, the system 
is strikingly different in other aspects and is described by 
a new non- Fermi liquid fixed point. A scal ing de scription 
of this fixed point was developed in Refs. 13611371 . 

While we largely follow the earlier work, we find that it 
is convenient to consider a slight reformulation in which 
the only uncontrolled approximation is to assume that 
the gauge field dynamics can be described correctly by 
the RPA approximation, retaining only terms quadratic 
in the gauge field with a singular quadratic kernel. A 
virtue of this approach is that one thereby obtains a the- 
ory which has an N flavor extension which is soluble at 
N = oo. Moreover, a controlled and systematic pertur- 
bation expansion in powers of 1/N can be implemented, 
which allows one to compute physical properties in terms 
of non-universal "bare" parameters such as the shapes of 
the Fermi surfaces. We employ this approach to calculate 
the leading 1/N behavior for both the boson correlator, 
Gfc(r), and the density-density correlator, D b (r), in the 
D-wave Bose Liquid. 

In the next Sec. IVIl we will consider an effective field 
theory approach, which allows us to check for the stabil- 
ity of the DBL phase that is present at large N. Specifi- 
cally, we study the effects of residual short range attrac- 
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tive interactions between the two fermion flavors. With 
vanishing D-eccentricity when the Fermi surfaces for both 
fermion species become the same, the s-wavc Cooper 
channel is "nested" , and a possible instability towards 
a paired BCS state seems likely. The resulting state is a 
bosonic condensate, since (&t) = (d\4) 0. However, 
a non-vanishing D-eccentricity of the Fermi surfaces de- 
stroys the BCS nesting and the possible pairing insta- 
bility. Indeed, we will find that the DBL with a large 
D-eccentricity can exist as a stable phase. (For smaller 
D-eccentricity, an instability towards a finite momentum 
Bose condensate or an incommensurate charge density 
wave co-existing with the DBL is a possibility.) As is 
discussed in Sec. IIV1 this is nicely in line with the prop- 
erties of the associated Gutzwiller wavefunctions: The 
(det) 2 wavefunction that obtains in the limit of zero D- 
eccentricity appears to have off-diagonal long range or- 
der, whereas the general (det) x (det)^ wavefunction ex- 
hibits power law correlations consistent with the DBL, 
as extracted from the gauge theory. 

A. Formulation 

Being interested in the low energy properties, it is le- 
gitimate to focus on the fermions living near the Fermi 
surfaces, just as in Fermi liquid theory. Moreover, the 
longitudinal Coulomb interactions mediated by the gauge 
field are readily screened out within an RPA treatment, 
generating short-ranged screened density-density inter- 
actions. But the transverse fluctuations of the gauge 
field - the photon - are incompletely screened by the 
fermion particle- hole excitations. It is thus adequate to 
just retain the transverse components of the vector po- 
tential. Working in the Coulomb gauge with V • a = 0, 
the vector potential reduces to a scalar, e.g., a(k) = 
o(k) sign(fc J/ )(fc J/ x - fc x y)/|k|. 

Moreover, for each fermion species, it will be sufficient 
to focus on a pair of Fermi surface patches with normals 
parallel to some axis, say x, as shown in Fig.[l0l It is the 
a x component of the vector potential which is minimally 
coupled to these fermions. The important wave vectors of 
the gauge field a x (k) that are strongly Landau damped 
by the fermions in these patches satisfy k x <C k v (see 
Fig. [TO]) , and in this region of momentum space a x (k) ss 
a(k). Conversely, the modes of the gauge field a(k) with 
k x <C k y feed back and scatter the patch fermions. Thus, 
we can focus on fermion fields c?i(k), (fe(k) and the gauge 
field a(k) that are confined to their respective patches in 
momentum space, as shown schematically in Fig. 1101 

For concreteness, we assume that for each fermion 
species there are two relevant patches on the opposite 
sides of the Fermi surface that are labelled s = R/ L = 
+/— for the right/left patch as indicated in Fig. [10] We 
then decompose the Fermion fields into right and left 
movers by writing, 

d Q (r)~£e iak '=- r d Qfl (r) . (51) 



Here, kp a denote the locations on the two Fermi surfaces 
which have normals aligned along the x-axis, and the 
fields d as (r) are assumed to be slowly varying on the 
scale of kp 1 . 

The full low energy action consists of three terms de- 
scribing the dynamics of the gauge field, of the fermions, 
and of their interactions: 

S<°> - J d 2 r dr [C d + 4°> + C mt ] ■ (52) 

The "patch" Lagrangian density for the fermions is sim- 
ply, 

C d = ^2 ^2 dl s [d T + sv a (-id x ) + -^(-id y ) 2 ]d as , 

a s— ± 

(53) 

which is characterized by the local Fermi velocities, v a , 
(a = 1,2), and the local Fermi surface curvatures, c a . 
The curvature is crucial in the fermion-gauge problem 
since there is an effective decoupling of the fermions out- 
side the patch from the gauge field, a = a x , which is 
parallel to the patch normals. Here s = ±1 specifies the 
sign of the group velocity while v will always denote the 
absolute value. 

The dynamics for the transverse gauge field describes 
the free photon, 

C^ = \[{d T af + K{Vaf], (54) 

whereas the gauge-fcrmion interactions are of the usual 
minimal coupling form: 

C int = ^2 g a sv a d) as a d as . (55) 

a, s 

Here, g\ = —g<i = g with g = 1, but it will be convenient 
for bookkeeping purposes to retain g as a dimensionlcss 
parameter in the theory. We emphasize that both the 
fermion fields and the gauge field, which enter the above 
action in real space, are slowly varying fields correspond- 
ing to momentum modes inside their respective patches; 
in what follows, we always specify the wave vectors of the 
slow fields relative to their respective patch centers. 

B. RPA approximation for the gauge field 
dynamics 

The above theory is strongly coupled and cannot be 
treated pcrturbatively in g. As such, it is necessary to re- 
sort to approximations to make any meaningful progress. 
Here we make our central approximation, namely re- 
placing the dynamics of the gauge field by a Landau 
damped form which results if one sums the RPA bub- 
ble of particle-hole excitations: S a — * S a with 

Sa = \J (fM |a(q,c)| 2 . (56) 
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FIG. 10: Relevant momentum patches of the strongly coupled 
fermion - transverse gauge field system. For simplicity, here 
and throughout the text, we take the Fermi surface patches 
with normals in the x direction; the required properties are 
the Fermi velocities and surface curvatures; the latter are cru- 
cial since differently oriented patch systems effectively decou- 
ple at low energies. 



The momentum/frequency integral is understood as 
/ = J d 2 qdui / (2tt) 3 , with the momentum restricted to 
a region satisfying \q x \ < A x ,\q y \ < A y with A x <C A y . 

The magnitude of the Landau damping coefficient, T, 
will be determined by the low-energy fermions near the 
specified Fermi surface patches. Within the RPA approx- 
imation one obtains, 



RPA 
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2tt 



1 

c-2 



(57) 



but we will retain T as an independent parameter since 
beyond this point we will not be using RPA in any event. 
It will sometimes prove convenient to consider the limit 
in which the dimensionless gauge charge, 



rv 



(58) 



is small. Here c is a characteristic Fermi surface cur- 
vature. Indeed, within the iV-flavor extension discussed 
below, while 1/iV will be the small parameter that leads 
to a controlled analysis, already at order 1/N there are 
some physical properties which we can only compute an- 
alytically when e 2 is simultaneously taken as a small pa- 
rameter. 

The parameter x that enters in S a is a "stiffness" which 
has to be positive for stability. Starting from the lattice 
gauge theory formulation, \ gives a measure of the en- 
ergetic cost of setting up static internal fluxes. One can 
very crudely estimate x using the expression for the dia- 
magnetic response of free fermions, 
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(59) 



where rrid is some effective mass. We reiterate that this 
stiffness comes from some energetics that is assumed 
to stabilize the Fermi surfaces of the d fermions in the 
first place. Therefore, from the low-energy perspective, 
X should be viewed as a phenomcnological parameter 
that encodes some high-energy physics. Furthermore, 



the stiffness x 1S independent of the patch orientation 
in momentum space for the transverse gauge field a(k), 
since gauge invariance requires that the field energy be 
proportional to xB 2 = x(^ a y — 9 y a x ) 2 . In any event, 
the magnitude of x will pl a y a minor role below. Indeed, 
for correlation functions with power law decay as we find 
below, changing % will only modify the amplitude but 
not the form of the correlators. 

Before proceeding with an analysis of the full theory, 
Sd + S a + Si n t, it is important to note that within the 
above approximation the Landau damped dynamics of 
the gauge field described by the action S a is both singular 
and harmonic. We have ignored possible terms in the 
action which involve higher powers of the gauge field. 
As such, at this stage we could integrate out the gauge 
field exactly generating a non-local and retarded four- 
fermion interaction. It is the absence of non-linearities 
for the gauge field dynamics which allows us to introduce 
a large- N generalization which is soluble at N = oo and 
which facilitates a formal and systematic 1/N expansion, 
as we now describe. 



C. N flavor extension 

Purely for technical reasons, then, we generalize the 
theory to N flavors of fermions at the two patches on 
each Fermi surface, d as — > d ias . Here i — 1, . . . , N labels 
the flavor index, while s — ± labels the patch location 
and a = 1, 2 the two Fermion species. The action for the 
fermions is simply generalized as, 



isv a d x 



VgCg n2\j 

Oy)(li as 



(60) 

The real scalar gauge field a(r) is likewise generalized, 
but now as an N x N Hcrmitian matrix of complex fields, 
cjjj(r) = a^j(r). When expressed in momentum space, 
each component of the matrix, a.y (q) , has dynamics with 
a Landau damped form, 

Si N} = lf:f (r^\+x\ qy \ 2 )K(^)\ 2 - (61) 
ij=l ->iu \ \ ( lv\ / 

Since ajj(q, lo) = a^-(— q, — u>), this can also be re-written 
as a trace over the flavor indices of the matrix *a*(q, oj): 



w)V(-q, -u))] . 

(62) 



Finally, the fcrmion-gauge field interaction is taken as, 

JV 



S 



' ^ Vn a .4- 



(63) 



Notice that the magnitude of the interaction strength 
(the charge) has been taken to scale as 1/yN. At N — 
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1 the theory reduces to that discussed in the previous 
subsection, with the assumption of a harmonic Landau 
damped form for the gauge field dynamics as the only 
approximation. 

It is worth pointing out that the full action has an 
enlarged global symmetry being invariant under, 

dias * ^ 6 ^ diets : &ij ^ ^ ^ ®ij ; (^^) 

for arbitrary phases Oi, i = 1, . . . , N, and (j) a , a = 1,2. 
Thus, the densities ^ Q d\ a d ia are independently con- 
served for any i, and also J2i dt a di a are independently 
conserved for any a. Note also that the dynamics of each 
of the A" flavor fcrmions dn is identical to one another, 
and similarly for all the di2 fermions. 

Before solving the full TV-flavor theory, S [ d N) + si N) + 

Sfat j exactly at N = oo, we define large N general- 
izations of the various correlation functions. As in the 
gauge mean field theory, we anticipate that the exact 
fermion Green's function for (any) one of the TV fermion 
flavors, which we denote G^\r), will be dominated at 
large distances by the patch fermions with normals paral- 
lel or antiparallel to the observation direction f . We can 
thus expand G^\r) for large |r| in terms of the "patch" 
fermion Green's functions as, 

^'(r.rl^^TGWj,.,,) , (65) 

S 

with the definition, 

GW(r-r',r-r') = (4 s (r, r)d ias (r', r')) . (66) 

Note that there is no summation over i; throughout, any 
summation over indices will be shown explicitly. 



There are now A" flavors of bosons with creation oper- 
ators, 

6t(r)=4(r)4(r), (67) 

with all N boson flavors having identical dynamics. The 
large N generalization of the boson correlator is then 
simply, 

Gg(r-r',r-r') = (6l(r 1 r)6 J -(r',r')) . (68) 

Due to the global symmetries, the correlators which 
arc off-diagonal in the flavor indices vanish, leaving 
a uniquely defined boson Green's function, = 

5ijGi N \ As desired, G^ reduces at A^ = 1 to the boson 
correlator defined in Sec. IIII Bl 

The boson density-density correlator, , is readily 
defined in terms of (any) of the boson densities pa as in 
Eq. (|34p . Similarly, we can define the fermion density- 
density correlators for each species, D^, in terms of any 
one of the fermion densities, pi a . As before, we assume 
that the boson density-density correlator can be approxi- 
mated as an average over a of the fermion density-density 

correlators, ~ \ ^ Q ■ 

In order to extract the above correlation functions, 
it is convenient to define patch fermion particle-hole 
and particle-particle bubbles, which we denote as IT^ 

and Tlp^\ respectively. The exact bubble correlators 
can be formally expressed in terms of the exact patch 
fermion Green's functions and the exact particle-hole and 
particle-particle vertices, 



J q'o;' 

I 



Here V^ h and Vp P ' denote the fully renormalized ver- 
tices, and a, s and a' , s' label the patch locations of the 
two fermion lines that come out of the vertex. The to- 
tal momentum and frequency running through the bub- 
ble is q, u), and is divided between the two fermions as 
shown above. At N = oo there are no vertex correc- 
tions, = 1 and Vp£°' = 1, and each bubble is just a 
convolution of the two Green's functions. 

The boson correlator at large distances and long times 
can be readily expressed in terms of the bubbles, 

G { b N) (sk Fl + s'k F2 + q, U ) » n^ ls;2s , (q, U ) , (71) 



and similarly for the fermion density-density correlator, 
D W {skFa s > kFa + q , u) « n W qs;qs , (q, W ) . (72) 

D. Solution at TV = oo 

At A" = oo the fermion Green's function G^ (q, uj) 
can be readily obtained by summing the nested rainbow 
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diagrams to give 



i 



with 



iuj (1 + |n a /«|V3) - SVa q x - *f-q> y 

(73) 



2^ v / 3x 2/3 r 1 / 3 



(74) 



Because of the gauge interactions, each fermion 
species obtains an anomalous self-energy of the form 
— iuj^l/uj^^ which dominates over the bare term — iuj 
on energy scales below ft. Physically, the fermions cease 
to propagate as free particles due to the strong scattering 
from the dynamical gauge field, becoming "incoherent" . 

We first consider the long distance and time behavior 
of G^\r,r). As in the mean field Sec. IIII1 the Green's 
function is dominated at large distances by the patches 
of the Fermi surface with normals along ±f . For the 
equal-time Green's function, we obtain, 



G^M) 



30F 



1/2 



1 



CW, (75) 



where the velocity v a and curvature c a characterize the 
d a Fermi surface patches with normals along ±f, and 
Gf/(r) is given in Eq. ([28]). On the other hand, the 
long time dependence of the local fermion Green's func- 
tion is not affected by the anomalous self-energy and is 
determined solely by the density of states vq at the Fermi 
energy, 



T 



(76) 



Turning to the physical correlators, at TV = oo there 
are no vertex corrections to the fermion bubble, so we 
have simply, 



4r ) w=-i4r ) wi s 



(77) 



which gives for the TV = oo equal-time density correlator, 

OL 1 1 

where D^ IF (r) is given in Eq. (j3~7| . Similarly, in the 
absence of vertex corrections at TV = oo in the particle- 
particle bubble, one has for the equal-time boson Green's 
function, 



16 V ^1^2 



1/2 



(79) 
(80) 



where G^ IF (r) is given in Eq. (|29|) . However, the time 
dependence of the local boson Green's function G(,(0,t) 
is unchanged from the mean field Eq. (|3"Tj) . Finally, the 



box correlator is given by Eq. (|47|) with appropriate TV = 
oo fermion Green's functions and decays with a power 
law envelope of —x~ s . 

Notice that both the boson and the boson density- 
density correlators at TV = oo fall off in space as |r|~ 4 , 
which is faster than their mean field counterparts. In- 
jecting a boson corresponds to creating two fermions. In 
the mean field, these fermions propagate independently 
and as free particles. On the other hand, in the present 
TV = oo theory that sums the nested rainbow diagrams, 
both fermions propagate "incoherently" due to their scat- 
tering from the gauge field, which effectively reduces the 
ability of the created boson to propagate, thereby leading 
to a faster decay of the boson correlator. At this level 
of approximation, the two injected fermions, while scat- 
tered by the gauge fluctuations, do not scatter off one 
another. The effects of interactions between the two in- 
jected fermions appear at order 1/TV, entering as vertex 
corrections. Since the two fermion species have opposite 
gauge charge, one expects that they will attract one an- 
other. As we shall see, though, they do not form a com- 
posite boson (a bound state, for example). But their mo- 
tion becomes strongly correlated allowing them to prop- 
agate more effectively when close by spatially. The net 
result, as we find below, is that the boson can propagate 
more effectively when the motion of the fermion pair is 
correlated, leading to a slower decay of the correlator 
compared to the TV = oo result. Thus, the effect on the 
boson dynamics due to the decoherence experienced by 
the fermions from scattering off the gauge fluctuations 
is compensated, perhaps only partially, by the attractive 
interaction between the pair of fermions. 



E. Fermion Green's function at order 1/TV 

The 1/TV contributions to the fermion Green's function 
are obtained by considering the nested rainbow diagrams 
and replacing one gauge field propagator G a (q,u>) = 
(r\id\/\q y \ + x\ly\ 2 ) _1 w ith the bubble correction SG a ~ 
jfG^Il 1 ^^. Formally, we can first calculate the rain- 
bows using G a + SG a everywhere and then extract the 
1/TV piece. As we now argue, the functional form of the 
fermion Green's function in Eq. (|73| remains unchanged 
except for finite 1/TV shifts in the parameters. Indeed, 
evaluating the TV = oo particle/hole bubble at small 



wave vector and frequency gives n^ 3 '' (q, uj 



\Qv\ 



which only leads to a finite 1/TV shift in T. This is 
not surprising, since the bare singular gauge propaga- 
tor postulated at the outset in our theory is motivated 
by the free fermion particle-hole bubble, and one can ver- 
ify that the bubble remains unchanged also in the pres- 
ence of an arbitrary fermion self-energy that depends on 
the frequency only. Then, using this renormalized gauge 
propagator to evaluate the rainbow diagram gives a con- 
tribution to the fermion Green's function of the form, 



2/3 [G { r } } 2 - This will at most give finite 
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1 /N shifts to il a , v a , c a . Next we consider the vertices 
at 1/N, which undergo more dramatic modifications. 



F. Vertices at order 1/N 



Single photon contribution 



Consider first dressing the bare vertices with a single 
gauge propagator. We denote the value of the corre- 
sponding Feynman diagrams as V^ a ' and Vpp , where 
(la) stands for one (Landau damped) photon exchanged. 
These diagrams are the first ones in Figs. [TT1 and [T2"l As 
we will discuss below, the full structure at 1/N is more 
involved, but the calculation of ]/( la ) already captures 
the main qualitative physics. 

For simplicity, we first focus on the situation with zero 
external y momentum and frequency, defining, 

v (ia) = ^(cbwjq'.o/Jl^^^^o , (81) 

where q, q', to, to' denote the momentum/frequencies in 
and out of the two legs and q~ = q x — q' x is the momen- 
tum running through the vertex. For notational simplic- 
ity we have suppressed the corresponding flavor labels, 
a, s; a', s' . Evaluating the Feynman diagrams gives, 



V 



(la) 



V 



(la) 



pp; <ix 



= jrSign(g a g al )S s - s 'X ph ln(A x /\q x \) , (82) 
= siga(g a g a/ ) 6 s - sl \ppln(A x /\q~\) (83) 



revealing a logarithmic divergence cut off by the total 
x-momentum running through the vertex. Explicit ex- 
pressions for the dimensionless non-negative coefficients 
X p h and X pp are given below. 

Notice that both the particle-hole and the particle- 
particle vertex corrections vanish unless the two patches 
have opposite group velocities. Moreover, the sign of 
the correction is given by the sign of the product of the 
two charges and differs for the particle-hole and particle- 
particle channels. The origin of these signs is essentially 
Ampere's law familiar from electrodynamics. Two mov- 
ing charges producing parallel charge currents experience 
an attractive interaction, bringing them closer together 
and enhancing their interaction strength, which is encap- 
sulated by the magnitude of the interaction vertex. For 
example, a particle and a hole of the same fermion species 



(a = a') residing on the opposite Fermi surface patches 
(s = —s') have an enhanced interaction strength, which 

is encoded in > 0. On the other hand, the in- 

teraction vertex is suppressed when the charge currents 
are anti-parallel. For example, a pair of particles with 
the same charge (g = g') experience a suppressed vertex 
interaction, Vpp < 0, when their group velocities are 
antiparallel. 

The absence of singular vertex corrections when the 
two patch fermions have parallel group velocities, s = s' , 
is, we believe, not just a peculiarity of this lowest or- 
der contribution, but will be valid generally for the exact 
vertices at arbitrary N. This expectation is based on the 
specific form of the propagators: the fermion propagator 
has simple poles in the complex q x plane while the gauge 
field propagator is independent of q x , the latter a reflec- 
tion of the transverse nature of the gauge field. For any 
given Feynman diagram, it is the sign of the internal fre- 
quency variables running around the loops together with 
s, s' which will determine whether the poles in the inte- 
grand will be in the upper or the lower complex q x planes. 
For vertices with s = s' all of the fermion lines internal 
to any diagram will have the same group velocity, and 
with energy (frequency) conservation one expects that 
at least one of the q x variables will only have poles in the 
upper (or lower) half-plane, so that the contour can be 
appropriately deformed to show that the diagram van- 
ishes. Thus, we henceforth restrict attention exclusively 
to vertices with s — —s', and to ease the notation will 
drop the explicit s, s' indices in the following. 

The dimensionless magnitudes of the vertex enhance- 
ments in the particle-hole and the particle-particle chan- 
nel, Xph and Xpp respectively, are given by, 



^Kph} , (84) 



^ph £[Cpfi<] i ^pp — ^[Cpp] • 



£[Cpp] ; (85) 



(86) 



In the last line, we specialized to the N = oo expressions 
for and fl' (Eq. l74"|) . which is valid for extracting the or- 
der 1/N vertex corrections. Here, £[£] is a dimensionless 
function of a dimensionless argument given by 



J 



3\/3 r°° tdt 

~2^J Q (l + t 3 )(l + C 2 t 4 ) 



3ttC 5 + 8tt(1 - C 4 )/\/3 - 3V2V(C 1/2 - C' 2 ) - 6C 2 In C 
8^(1 + C 6 )/V3 



(87) 



and has been normalized so that £[0] = 1. This func- tion is monotonically decreasing with increasing argu- 
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FIG. 11: Ladder diagrams that contribute to V p h at 1/N. 





FIG. 12: Crossed diagrams that contribute to V pp at 1/N. 



ment varying as £ [Q sa ^1 for ( 
C p h and C PP arc given by, 



Cph 

Cpp 

Cph 



r 2 / 3 (c 



2x2/3 + ^ 

r 2 / 3 |c- c '| 



oo. The parameters 

(88) 
(89) 
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V3tt 
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r|c-c'| . (90) 



In the last line, we again used the N — oo values for fl, O'. 
Notice that the vertex interaction strength between two 
fermions on opposing patches depends on their particular 
Fermi surface curvatures c, d . Moreover, in the particle- 
particle channel the interaction strength is maximal when 
the curvatures are equal. This reflects a particle-particle 
nesting, which is responsible for the BCS pairing insta- 
bility in the Fermi liquid context. 



2. Photon ladders 

Having evaluated the one-photon corrections y( la ) to 
the bare vertices, we now discuss the full set of diagrams 
that appear at order 1/N in each case. First consider 
the fcrmion particle-hole vertex which enters into the 



density-density correlator: V f 



(N) 



V. 



ph; a,a " 



The 1/N 



diagrams are shown in Fig. QT] and contain an arbitrary 
number of non-crossing photon lines,— beginning with 
the one-photon diagram V^ a \ Eq. (|82[) . We shall de- 
note the contribution from the diagram with n photon 
lines as V^ a) . The magnitude of the singular contribu- 
tions from each diagram in the ladder sum will depend 
on a dimcnsionlcss gauge charge, 



.9 2 /rc Q 



(91) 



Above we found V. 



(la) 



ph 



(l/N)£[V3n/el)Hl/q- 



(1/N)A^ \n(l/q x ). The singular contributions from the 
diagram with two photon lines can also be extracted from 
the Feynman diagram, V 



A£' ln^(l/g~)], with Af'jA^' universal dimensionless 
functions of e 2 . From this we infer the general structure 
for the diagram with n photon lines, 



(2a) 
ph 



= (l/N)[A^Hl/ q -) 
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(92) 



This form suggests that the full set of 1/N ladder dia- 
grams can be exponentiated, 



1 

N 



exp[j a \n(l/q x ; 



with 7 Q = 2™ i ^4-1 • The validity of this exponentia- 
tion can be justified by formally summing the ladder se- 
ries to obtain an integral equation for the full vertex. The 
integral equation has a singular kernel and its solution is 
expected to Imvv the power-law form with some exponent 
7 Q which will depend on the dimensionless charge e 2 . 

Unfortunately, we have been unable to solve the inte- 
gral equation to obtain j a . However, we observe that for 
small charge e 2 , a[ 1] = 3e 2 /(87r) + 0(e%). This gives 
the leading behavior even though A^ is not analytic in 
e 2 , as can be seen from Eq. (|87[) ; more generally, one 



1 

N 



\9x - q x \ 



(93) 



finds A 



(n) 



0(el 



Thus, assuming small parameter 



e 2 , the leading contribution to 7 Q comes just from the 
single photon diagram, 



7a 



-el + 0(e 4 a ,l/N) 



(94) 



For e 2 of order one, the exponent j a , while obtained 
by summing diagrams at order 1/N, is itself of order 
one. It is the amplitude of the vertex in Eq. |93|) which 
is of order 1/N. This reflects the particular structure of 
the large N theory, which, constructed to sum the nested 
rainbows for the fermion self energy at leading order, also 
exponentiates the ladder sum for the particle-hole vertex 
at order 1/N. Notice that j a monotonically decreases 
as the curvature of the Fermi surface is taken large, due 
to the increasingly restrictive phase space requirements 
on the particle- hole vertex in this limit. We expect that 
this trend is valid more generally, for larger e 2 and for 
arbitrary N. 

Summarizing the discussion of the particlc-hole vertex, 
we can also write a scaling form for more general external 
momenta and frequencies, 



V p W~\ qx - q > x \-^V Pi 
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(95) 



We next consider the vertex which enters in the boson- 



boson correlator. V, 



(AO 



V, 



(AO 



The full set of di- 



agrams that contribute to 



pp; 1,2- 

the particle-particle vertex 
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are shown in Fig. [12] and contain an arbitrary number 
of maximally-crossed photon finest 2 - The first term is 
the one-photon diagram that we evaluated in Eq. ([83]) . 
V£$, 2 = (l/N)S[C PP ]Hl/q-) = (l/JV)5iln(l/g-) 
(plus non-singular contributions). An evaluation of the 
two-photon diagram gives the same form, 2 = 

(\/N)B2 ln(<7~). In contrast to the particle-hole case, 
here the two-photon diagram does not contain a log- 
squared term, which indicates that the crossed ladder 
series cannot be exponentiated. While we cannot evalu- 
ate the higher order diagrams, we expect that each term 
will have a log contribution with some dimcnsionlcss am- 
plitudes B n . If this is the case, the full structure of the 
particle-particle vertex up to order 1 /N is, 



G. Boson and density correlators at finite N 

Since there are no singular 1/N corrections to the 
fermion Green's functions, to extract the boson and den- 
sity correlators up to order 1/N one can use the N = oo 
fermion Green's functions when evaluating the bubble di- 
agram. Using Eq. ([75)1 and the scaling form Eq. ([9"9"| for 
the Vb vertex enables the particle-particle bubble to be 
expressed as, 



|-(2+»7) 
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where P is an appropriate scaling function. One can then 
scale out the q x dependence under the integrals to obtain, 



(N) 



= 1 



rji ln(A x /q x 
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(96) 



n^)(q, w )~| & | 1 -"n p 



q x ' q x 



(100) 



with rji = Yl'n—i B n - To check whether the latter series 
can be defined in principle, we again take the limit of 
small dimcnsionlcss charges e 2 . Consider the least favor- 
able case with matched Fermi surface curvatures, c\ = C2, 
so e? = en = i 



We thereby arrive at the desired expression for the boson 
correlator. 



g 2 /(Tc). We find that the first term is (k Fl - kp 2 + q, w) ~ l&I^G I ^, 



.2/3 



(101) 



independent of e 2 , B\ = 1. However, evaluating the two- 
photon diagram gives B2 ~ (e 2 ) 2 ln 3 (l/e 2 ) for small e 2 . 
The diagrams are expected to be smaller when C\ 7^ c 2 
(for example, even the first term B\ = £[Cpp] can be made 
small when g 2 <C r|ci — C2I). We thus conjecture that 
the series converges for small enough e 2 . 

Provided it is legitimate to exchange the order of limits 
N — > 00 and q x — ► in Eq. ([96]) . the boson vertex can be 
written in a power law form, 



KW(q,u,;q>' 



1 1y =q' y =u=w 
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with 



r) = m/N + 0(l/N 2 ) . 



(97) 



(98) 



We will discuss the legitimacy of this exponentiation pro- 
cedure in Sec. IVII For more general external momenta 
and frequencies the full boson vertex will satisfy a scaling 
form, 



V/, (JV) (q, W ;q>') 
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(99) 

where we have used the fact that the vertex depends 
only on the total x-momentum running through it, q~ = 
q x —q' x . The scaling function will vary as, Vb(X, 0, 0, 0) ~ 
X~ v when any of its arguments X — > 00. For example, 
one has, V b - \q y \- 2r > ~ |w|~ 2? '/ 3 . 

As will be detailed shortly, the singular vertices 
Eq. (jnS]) and Eq. ([9^|) will lead to anomalous decay ex- 
ponents for the boson density-density correlator and the 
boson Green's function. 



with exponent rj given in Eq. |98|) . While this result was 
obtained systematically to order 1/N, we expect it will 
remain valid more generally, and in particular will involve 
the same scaling combinations q 2 /q x , Lu 2 ^ 3 /q x . The scal- 
ing function G and the dependence of the exponent 77 on 
the bare parameters such as the Fermi surface curvatures 
will presumably vary with TV. Near the other momentum 
kp x +ki? 2 , the boson correlator will satisfy a similar scal- 
ing form except with rj = 0. As argued before, we expect 
that the anomalous exponent at momentum kp^ + k^ 9 
will vanish in general, not just to leading order in 1/N. 

The large distance behavior of the equal-time bo- 
son Green's function in the DBL follows directly from 
Eq. (| 1 1 [) and was already listed in the Introduction, 
Eq. ([5]). We remark that simple "power counting" using 
Eq. (jlOip works here because the relevant frequencies and 
momenta arc indeed in the "scaling regime" ; however, it 
should be kept in mind that we are dealing with line 
singularities and should exercise more caution with such 
arguments in general. Once again, both the Fermi wave 
vectors k^ 17 ki? 2 and the scaling exponent rj will depend 
on the location and shape of the Fermi surface patches 
with normals along the observation direction f . At wave 
vector kpj — kp 2 , the Green's function decay becomes 
more slow due to the Ampercan attraction between the 
d\ and e?2 fcrmions. The fermions tend to move as pairs, 
but unlike a Cooper pair their motion is not phase coher- 
ent - the pairs are uncondensed. The boson correlator, 
while decaying more slowly with positive 77, does not ex- 
hibit ODLRO in the DBL phase. 

As we found to leading order in 1/N, Eq. ([9"5[) . 7] is 
expected to be largest when the two Fermi surface curva- 
tures are equal. In the DBL phase which has closed Fermi 
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surfaces, the characteristics of the d\ and di patches with 
normals along either diagonal x ± y are always equal due 
to the symmetry of the square lattice. Moreover, the 
decay of the boson Green's function Gfc(r) will be non- 
oscillatory along the diagonals, 



G b (r) 



1 



r II (x±y) , 



(102) 



a prediction that we find consistent with the properties 
of the (det) x (det)j / wavefunctions studied in Sec. IIVI (see 
Fig. [7]). For c\ — C2, a very rough estimate can be ex- 
tracted from the leading behavior, rj = (1/N) + 0(e 2 /N) 



by setting N 



1 giving rj ss 1 . This would shift the 



N = oo decay exponent back to its mean field value of 
l/|r| 3 . As noted earlier, the gauge fluctuations encapsu- 
lated already at N = oo strongly modify the motion of 
each fermion leading to a boson that moves less coher- 
ently decaying with a larger exponent than in the mean 
field theory. The Ampercan attraction between the pair 
of fermions which enters at order 1/N compensates this 
effect, leading to a slower decay of the boson correlator. 
For r\ = 1 the two competing effects exactly compensate 
one another. 

Next we consider the density-density correlator. Just 
as for the boson correlator, the density correlator will 
satisfy a scaling form: 



D { d N) (2k Fa +c ll iu)^\q x \ 1 -^D 
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(103) 

with an anomalous exponent that depends on the Fermi 
surface patches, see Eq. (|94j) . Back in real space, the 
dominant density correlator oscillates with wave vector 
2kp a , with an envelope decaying as a power law with the 
anomalous exponent, 



cos[2k Fa -r-37r/2] 



ii ■ (104) 



the Fermi surface. Similar behavior was found for the 
density-density correlator. Provided this qualitative be- 
havior persists down to the physically relevant case of 
N = 1, we conclude that the DBL phase is described by 
a manifold of scale invariant theories rather than an iso- 
lated fixed point. Within the 1/N expansion, the power 
law form obtained for the boson correlator relied on the 
exponentiation of the logarithmic behavior on momen- 
tum and frequency of the leading 1/7V correction. In the 
next section, we check the legitimacy of this procedure 
by revisiting the renormalization group approach devel- 
oped earlier to describe the low energy physics of a sea 
of fermions coupled to a U(l) gauge field. 

We conclude this section by pointing out that through- 
out we worked on the assumption of having both d\ and 
e?2 parallel patches present as in Fig. [TO] In the case 
with open Fermi surfaces shown in the right panel of 
Fig. [TJ we encounter situations when for a given obser- 
vation direction f only one Fermi surface (or even none 
at all) has patches with normals in this direction. In 
such cases, the boson correlator will decay exponentially; 
there is no inter-species gauge interactions (one species 
is simply absent), but the fermions that are present are 
still strongly affected by the gauge field, and in particu- 
lar the preceding analysis of their density correlations 
remains the same. In the DLBL phase, there are no 
parallel patches on the entire two Fermi surfaces, so the 
d\ and d2 fermions effectively decouple from each other. 
It is also interesting to remark that in the limit of ex- 
treme D-eccentricity when the Fermi surfaces are com- 
pletely flat, we expect that the relevant gauge field is 
very strongly damped. This can be seen, e.g., from the 
naive RPA approximation for the Landau damping coef- 
ficient r, Eq. (|57|). in the limit of vanishing curvatures. 
In this limit, the energy scale f£, Eq. (|74|) . below which 
fermions become incoherent, goes to zero, which leads 
us to speculate that perhaps in this case the d\ and c?2 
fermions behave as essentially free. 



We have also indicated that the zero momentum com- 
ponent is not modified relative to the N = oo behavior. 
Again, both hp a and 7 Q are particular to the Fermi sur- 
face patches with normals parallel to ±f . In contrast 
to the particle-particle channel, the particle-hole channel 
is not nested so that we suspect the exponents 7 Q will 
be smaller than rj for the physically relevant case with 
N = 1 . 

In sum, at leading order in our systematic 1/N ex- 
pansion, the boson correlator was found to exhibit power 
law decay with an exponent that depends on the bare 
Fermi surface curvatures, varying continuously around 



VI. EFFECTIVE FIELD THEORY FOR THE 
DBL: CRITIC ALITY AND STABILITY 

Polchinski 3 ^ and other a 37 ' 38 argued that 2d fermions 
interacting with a U(l) gauge field can be fruitfully stud- 
ied within a renormalization group analysis of a partic- 
ular effective field theory. In the present case, the fixed 
point theory has an action given by, 



Sfixed point — Sd + S a + iSj, 



(105) 



with 



as ' 



s2 Q 



1/3 



v a c, 



+ sv a q x H — q d as (q,uj) , 



(106) 
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while the gauge field part S a is given in Eq. and the 
fermion-gauge coupling is specified in Eq. (|55[) . Within 
this fixed point ansatz, both the gauge field and the 
fermions have singular propagators, but their interac- 
tion is local. The advantage of this effective field theory 
is that it is possible to define a simple renormalization 
group transformation which leaves the Gaussian part of 
the theory invariant, and under which the interaction 
strength g is a marginal perturbation. If g is assumed to 
be small, it can be treated via a conventional perturba- 
tive RG approach as we describe below. One finds that 
the full effective field theory is invariant under the RG, 
i.e., it is at a fixed point. This strongly suggests that ex- 
ponentiating the logarithmic singularities present in the 
large N approach, as we did to order l/N above, is a 
correct procedure at all orders in l/N. If it is, then at 
order l/N 2 one must find (log) 2 singularities with partic- 
ular coefficients such that the exponentiation procedure 
to obtain a power law behavior is consistent. 

An unsettling drawback with this effective field the- 
ory approach is that it leaves unclear what constraints 
must be placed on additional interactions that can be 
added to the theory. At the very least these interactions 
must be gauge invariant and consistent with momentum 
conservation, but is it legitimate to require that the inter- 
actions be local in the fcrmion fields? For example, can 
one require that four-fermion interactions be local with 
non-singular coefficients, or will singular interactions nec- 
essarily be generated? Indeed, if one were to formally 
integrate out the gauge field with its singular propaga- 
tor, one would generate four-fermion interactions with 
a particular singular form. In what follows we will ig- 
nore these subtleties, exploring the possible perturbative 
instabilities driven by non-singular fermion interactions. 
Specifically, we will consider all non-singular quartic in- 
teractions involving four fermions living near the Fermi 
surfaces that are consistent with the relevant symmetries, 
most importantly momentum conservation. 

The fixed point action, as it stands, has 9 parameters, 

n a ,g. (107) 

But with an appropriate rescaling of the fermion fields 
and the gauge field together with the momenta and fre- 
quency, it is possible to set 5 of these parameters to unity, 
T = x = g = v a = l- The effective field theory is then 
specified by 4 dimcnsionlcss parameters, c Q = c a T/g 2 
and £l a = VL a x 2 T /{g 2 vl/). Evidently, this theory is not 
describing a fixed "point" per se, but constitutes a four- 
dimensional manifold of theories which are invariant un- 
der the RG. Establishing definitively that a particular 
bare (lattice) gauge theory Hamiltonian is attracted to 
this manifold is exceedingly difficult. Arguably, it is 
even harder to deduce where on this manifold the theory 
flows. The values of the dimcnsionlcss parameters which 
are obtained from the leading large N analysis can, nev- 
ertheless, be used as a rough guide in addressing both 
questions. 



The RG analysis of the fermion-gauge action, 
Sfixed point , proceeds as follows. At each stage, the 
fcrmion fields reside in the appropriate momentum space 
patches l^l < A, \q y \ < (A/c) 1 / 2 , where A is the 
shell width in the direction normal to the Fermi sur- 
face. The corresponding restriction on the frequency is 
\uj\ < (v A) 3 / 2 /^ 1 / 2 . The gauge fields reside in similar 
momentum-frequency regions but centered around zero 
momentum, and the overall setup is illustrated in Fig. 1101 
If the dimcnsionlcss parameters exhibited earlier are of 
order one (which is the case if the parameters are taken 
from the large- N analysis) , the corresponding regions are 
roughly similar for all fermion and gauge fields. Also, in 
practical RG calculations beyond the tree level it is con- 
venient to keep the cutoff only on the frequencies and 
perform unrestricted integrations over the momenta. 

We integrate out the high-energy fields from the shell 
between A and A/6 and then rescalc the momenta and 
frequencies in order to recover the initial cutoff: q x = 
q' x /b, q y = q'y/b 1 / 2 , and uj = lo' /b z / 2 . We also perform ap- 
propriate rescaling of the fields: d — b 2 d! , a — b 2 a' . Upon 
such tree-level scaling, the fixed-point action remains ex- 
actly as before, i.e., all couplings remain unchanged. It is 
also useful to bear in mind that even though we nominally 
restore all cutoffs as we proceed, in terms of the original 
momenta the patches become more and more elongated 
in the y direction in Fig. 1101 In particular, any overlap 
between two patches that correspond to nearby but non- 
parallel tangents y and y' goes to zero in the low-energy 
limit. 

To proceed with the analysis beyond the tree level, we 
will work perturbatively in the (dimensionless) fermion- 
gauge coupling, assuming that it is small. Quite gener- 
ally, since the momentum shell RG cannot produce terms 
that are singular at small frequencies and momenta, the 
couplings £l a and T will not renormalize at any order 
in the perturbation expansion. On the other hand, the 
parameters v a , c Q , \i an d 9, can potentially flow. How- 
ever, a direct examination shows that they do not renor- 
malize at the lowest (one-loop) order, and we strongly 
suspect that this remains true at all orders. The im- 
plication is that the fermion-gauge interaction is exactly 
marginal and that Sfixed point describes a manifold of RG 
fixed points parameterized by 4 dimcnsionlcss couplings. 
Although the fermion-gauge coupling need not be small, 
we expect that this RG-invariant manifold will extend 
outside the perturbatively accessible regime. This mani- 
fold describes the putative DBL phase. 

In order to establish the stability of the DBL, how- 
ever, we need to consider the effects of all symmetry 
allowed perturbations that can be added to Sfixed point • 
Stability requires that all such perturbations are irrele- 
vant under the RG. Restricting ourselves to local terms, 
we focus now on the four-fermion interactions which are 
most likely to destabilize the fixed point manifold. 
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FIG. 13: One-loop diagrams needed to calculate renormaliza- 
tion of four-fermion interactions. 

A. Short-range fermion interactions 

Let us consider a general quartic term, 

W di iSi dl 2S2 d a3S3 d a4S4 , (108) 

which can contain different fermion species labelled by 
a and involve different patches labelled by s. The 
fermion/patch labels are assumed to be such as to satisfy 
the momentum conservation. Also, we do not show ex- 
plicitly the (conserved) frequency and momentum of the 
fermion fields. The amplitude for this four-fermion inter- 
action, W, measures the strength of a direct scattering 
between two fermions. 

At the tree level in the RG, all such short-range inter- 
actions are irrelevant, 

W -► W = W/b , (109) 

scaling towards zero at low energies. Physically, the in- 
coherent motion of the fermions, manifest in their uj 2 ^ 3 
self energy, leads to a reduced phase space for direct two- 
body interactions relative to the case for free fermions. 

At one loop, a gauge propagator can connect any two 
legs of the quartic vertex which involve fermions located 
on patches with parallel or antiparallcl Fermi surface nor- 
mals. Two such physically distinct diagrams are shown in 
Fig. [131 Both processes correspond to an interaction be- 
tween the two fermions mediated by the Landau damped 
gauge field. The one-loop integrals needed to evaluate 
these diagrams are identical to those for the vertex cor- 
rections in subsection lV F 11 and give the following contri- 
butions to the four-fermion interaction before any rescal- 
ing: 

a) : SW = sign(g a g a ,)S s ,- s ,\ ph Wki(b) , (110) 

b) : SW = -sign(g a g a ,)5 s ,- s ,\ pp Whi(b) . (Ill) 

Here, the primed and unprimed parameters refer to the 
possibly distinct species of the two fermions connected 
by the gauge propagator. The anomalous exponents X p h 
and X pp are given in Eqs. ([53)) and (|53|) - these formu- 
las are general and do not assume any relations among 
the parameters, as appropriate in the present RG setup; 
sometimes, which will be indicated explicitly, we will use 
the N = oo expressions Eq. (f8l>)) to get a crude estimate 
of the numerical values. 



Similarly to the analysis in Sec. IV F 11 either diagram 
gives zero if the fermions involved have parallel group ve- 
locities, whereas a nontrivial contribution is found when 
the two fermion fields reside on the opposite patches with 
antiparallel group velocities. The sign is again deter- 
mined by Ampere's law: Two fermions emerging from a 
scattering process with parallel charge currents attract 
one another, spend more time close together, and conse- 
quently enhance the amplitude for the scattering strength 
W. Conversely, the interaction strength is suppressed 
when the pair of scattered fermions have anti-parallel 
charge currents. The difference in the absolute values 
of the contributions from the diagrams a) and b) is sim- 
ilar to the distinction between the particle-hole (CDW) 
and particle-particle (BCS) processes in a Fermi liquid. 
Thus, in the situation when c and c! have the same sign, 
the particle-hole process is not "nested" because of the 
curvature of the Fermi surface. On the other hand, the 
particle-particle process is somewhat better nested and 
becomes perfectly nested when c = c', which is the fa- 
miliar phase space reasoning in the theory of the BCS 
instability. This distinction between the processes can be 
clearly seen from the arguments C, v n and £ pp , Eqs. (J5SJ 
and (|89p . of the function £[(] in the two cases. 

With the above general expressions in hand, we now 
consider specific four-fermion interactions. Of interest is 
the one-loop contribution T oc g 2 to the eigenvalue of 
such perturbations at the fixed point, 

^ = (-l + T)W + 0(W 2 ) . (112) 

The 0(W 2 ) terms, which in general mix different four- 
fermion terms, can be deduced similarly to the RG treat- 
ment of quartic interactions around the Fermi liquid fixed 
poin t 43 i 44 and in fact have very similar structure. In the 
Fermi liquid case, the interactions are marginal and there 
is no term linear in W, so the 0(W 2 ) terms determine the 
physics. In the fermion-gauge case, on the other hand, 
the interactions can be either relevant when T > 1 or 
irrelevant when T < 1. For T < 1, the DBL phase is 
stable to weak perturbations, but sufficiently strong bare 
interactions might still drive the DBL through a quan- 
tum phase transition into another phase. In this case the 
0(W 2 ) terms could be helpful in discerning the nature 
of the new phase. 

As in Fermi liquid theory, the strongest constraint on 
the allowed quartic interactions comes from momentum 
conservation, since all four fermions have to live near the 
Fermi surfaces. In Fermi liquid theory there are only two 
allowed types of vertices, the forward scattering interac- 
tions which contribute to the Landau parameters and the 
Cooper vertices. We first consider interactions in which 
all four fermions are of the same species. A general for- 
ward scattering interaction then takes the form, 

<k'4Aa4cA' a , (H3) 

where k, k' are two wave vectors on the species a Fermi 
surface. Consider first the generic case with k ^ ±k'. 
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Since this interaction does not involve opposite patches 
on the Fermi surface, the one-loop contribution will van- 
ish. = 0. This is also the case when k = k'. The 
forward scattering vertex with k = — k' is part of the 
Cooper channel to which we now turn. 

For a single species of fermion, there exists only the 
odd angular momentum Cooper pairing channel ( "triplet 
pairing"), 

W£,k'4*<£kad-k'adk'a, (114) 

with W£ k , = — Wl k k , = — W k _ k ,. The one loop con- 
tribution is negative, T k k , = — (X pp + X pp ) for k ^ ±k', 
due to the Amperean repulsion between the two parti- 
cles (and between the two holes) with antiparallel group 
velocities. For the special case k = ±k' one has T k ±k = 
— 2(X pp — X ph ) since there is now also an Amperean at- 
traction between particles and holes on opposite patches 
of the Fermi surface. But since we expect X pp > X ph due 
to the nesting in the particle-particle channel, one has 
T* < in this case as well. Thus, the one-loop contri- 
bution makes this triplet pairing vertex more irrelevant. 
Physically the repulsive gauge interaction between two 
d a fermions is unfavorable for their pairing. 

We now turn to vertices which involve two d\ fermions 
and two d,2 fermions. For simplicity, we first focus on 
the situation with vanishing D-ecccntricity, so that the 
two Fermi surfaces coincide. Consider first the forward 
scattering interactions coupling the densities of the two 
fermion species, 

<;k; 2 4Ai4'2^' 2 . (115) 

The only non-vanishing one loop correction is when k = 
-k', and gives t£££ = 2(A k p - A k h ) > 0. If this is 
greater than one, it would signal a possible instability, 
but in any event we will see that this is smaller than 
the eigenvalue in the conventional BCS "singlet" pairing 
channel that we consider next. 

The singlet BCS interaction is of the form, 

W / k,k'4i d -k2 d -k'2rfk'i • (H6) 

With zero D-eccentricity as assumed here, the BCS in- 
teraction is nested. Moreover, since the two particles 
(or the two holes) here carry opposite gauge charge and 
therefore have parallel gauge currents, there will be an 
Amperean attraction giving a positive contribution to 
the Cooper vertex: T k k , = X pp + X pp , for k ^ ±k'. 
For a circular Fermi surface, symmetry dictates that X pp 
will be k- independent, but this is actually the case gen- 
erally for such one-loop contribution in the case with 
matched d\ and d% Fermi surfaces since the curvatures 
on the opposing patches are equal. Thus, as claimed ear- 
lier, the eigenvalue of the singlet BCS pairing interaction, 
T s = 2X pp . is larger than the corresponding eigenvalue 
of the forward scattering interaction coupling the den- 
sities of the two species on opposite sides of the Fermi 



surface, T s > T k ^ k . For the special case k = k', 
there will be an additional positive contribution to the 
eigenvalue of the singlet BCS pairing interaction due to 
the Amperean attraction between particles and holes, 
T k _ k = 2(A k p + Xp h ). But being a set of measure zero, 
this will not contribute to a BCS-type pairing instabil- 
ity. At the present it is not clear what state would be 
preferred by a large interaction of this form and what 
significance this might have. 

It is interesting to use the N = oo parameters to de- 
duce an approximate value for the one-loop contribution 
to T s . With zero D-eccentricity, the two Fermi sur- 
faces coincide, so the curvatures are equal, c\ = C2, and 
A k p = £ [0] = 1 for all k. This gives X s = 2 and implies 
that the BCS interaction is strongly relevant, 

dW s 

— = (-l + T s )W s = W s , (117) 

within this approximation. Since we expect some bare 
short-range attraction between the d\ and c?2 particles, 
the runaway flows will lead to BCS pairing. The com- 
posite boson = d\d\ will condense, resulting in a su- 
pcrfluid phase. Thus, based on this analysis, we suspect 
that the putative S-type Bose liquid phase accessed from 
the gauge theory with identical Fermi surfaces for d\ and 
the di will be gcncrically unstable towards superfluidity 

We now consider the situation with non-zero D- 
ccccntricity. In this case, d\ and di have different Fermi 
surfaces and the Cooper vertex will no longer be nested. 
This precludes a conventional weak coupling BCS insta- 
bility for the D-wave Bose liquid phase. The remaining 
channels that could potentially drive a weak coupling in- 
stability are the forward scattering interactions given in 
Eq. pi5p with the locations k and k' on the two Fermi 
surfaces tuned to have anti-parallel patch normals. 

For bosons on the square lattice, with increasing D- 
eccentricity the two closed Fermi surfaces will eventually 
open up, and above some critical eccentricity the require- 
ment of anti-parallel d\ and di Fermi surface normals 
will be no longer possible. In this large D-ecccntricity 
regime, when the two Fermi surfaces are so dissimilar, 
there are no such forward scattering interactions coupling 
the two species which can be enhanced by the gauge fluc- 
tuations. Sometimes one can find interactions that can 
be enhanced via the particle-hole attraction within the 
same species, e.g., d k i^-ki^k'2^k"2! which requires tun- 
ing both k' and k" to satisfy momentum conservation 
but this is not always possible. In any event, the corre- 
sponding one-loop T = X pfl is not likely to change the ir- 
relevance of this term from the tree level; e.g. the TV = oo 
estimate of A k /i; Eq. (|86p . is always smaller than 1. Re- 
garding other interactions, none of the intra-species four 
fermion terms have positive eigenvalues. We conclude 
that in this large D-eccentricity regime the D-wave Lo- 
cal Bose liquid (as specified in Sec. IIII C[) will exist as a 
stable phase. 

Returning to the case with smaller D-eccentricity, the 
stability of the DBL phase will depend on the eigenval- 
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ues T-^(n) of the forward scattering interaction that cou- 
ples a patch on one Fermi surface with normal n to a 
patch on the other Fermi surface with anti-parallel nor- 
mal. The particles on opposing patches experience a 
strong Amperean attraction, whereas a particle and hole 
experience a weaker repulsion, so that the eigenvalues 
T-f (n) = 2(X P p — X p h) will be positive for all orientations 
of n. However, the magnitude of T-^(n) depends on the 
Fermi surface curvatures and will vary with orientation. 
This can be made explicit by evaluating the one-loop 
contribution using the N = oo values of the various pa- 
rameters, thereby we obtain an approximate expression 
for the eigenvalue, 
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The first contribution is maximal and equal to 2 when 
the curvatures on the two Fermi surfaces coincide, which 
they will when n is along a diagonal of the square lattice. 
We have no such estimate of the second contribution since 
g, r, c a are independent in our N = oo theory; if we use 
in addition a crude RPA approximation Eq. (|57j). we get 
the total T f = 2- 2£[VS\ = 1-53 > 1. In this case, the 
D-wave Bosc liquid in this regime of D-eccentricity will 
presumably be unstable, driven by the forward scattering 
interactions with n = x ± y along the diagonals of the 
square lattice. 

The above numerical estimate is very uncontrolled and 
is given only to see possible trends. Still, let us specu- 
late what the resulting phase might if there is indeed 
such instability. Most naively, there are two guesses, ob- 
tained by a mean field decoupling of the four fermion 



interaction dl, <iki<iL„dk'2, where k and k' are the loca- 
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tions of the two patches with normals n and — n that 
give the largest We can assume a non- vanishing 

particle-particle or particle-hole condensate, (dj. dL, 2 ) or 

(<i kl dk'2), respectively. Both condensates carry non-zero 
momentum, k + k' for the particle-particle condensate, 
and k — k' for the particle-hole. But since the inter- 
action is repulsive in the particle-hole channel, such a 
condensate seems rather unlikely. Also, the combination 
e?{e?2 carries a non-zero gauge charge, so is unphysical 
and cannot serve as a legitimate order parameter [but 

the product (c?ki^k'2<^lk'2^-ki) is g au S e invariant, and 
if condensed would correspond to an energy density wave 
at momentum 2(k — k')]. 

If there is an instability, we think it is more likely to 
occur in the particle-particle channel. The order parame- 
ter (rfki^k'2) i s S au S e invariant and correspond physically 
to a finite momentum Bose condensate. The situation is 
analogous to the FFLO problem, and as there, provided 
the condensate is not too strong, it will only gap out 
parts of the Fermi surfaces. As such, this would still be 
a very unusual type of superfluid, which in a mean field 
description would have residual gapless fcrmionic exci- 
tations. Beyond the mean field, the gauge fluctuations 
would still scatter the gapless fermions, rendering them 



incoherent. This state would thus correspond to finite 
momentum Bose condensation co-existing with the DBL. 

In sum, whether or not the DBL for small D- 
eccentricity is unstable to such finite momentum pairing 
is a quantitative issue which will depend on the values 
of the anomalous dimensions. An analysis of the actual 
wavefunctions in Sec. IIVI reveals no such tendencies, so 
might be taken as an indication against the instability. 
Indeed, as we saw in Sec. IIVI the S-wave Bose liquid 
wavefunction does appear generically unstable towards a 
conventional superfluid, consistent with the expectation 
of a zero momentum BCS instability in the gauge theory. 
In the same spirit, we will see in the next section that 
for hard-core bosons at half-filling the DBL wavefunc- 
tion in the extreme D-eccentricity limit (corresponding 
to fermions which can only hop in one of the two di- 
rections on the square lattice) reveals a co-existence of 
a commensurate (ir, ir) CDW with a gapless DBL Fermi 
surface. This indicates that the Gutzwillcr wavefunc- 
tions for the DBL, at least in some special instances, can 
also reveal translational symmetry breaking instabilities, 
while we repeat that no such instabilities are observed 
for generic DBL wavefunctions. 

This concludes our discussion of the effective field the- 
ory description and possible instabilities of the DBL 
phase. As we have seen, the main potential instability 
involves pairing d\ and d,2 fermions moving with opposite 
group velocities, but such BCS channel is suppressed for 
mismatched Fermi surfaces and is completely eliminated 
in the DLBL regime, which appears to be particularly 
stable. 



VII. RING HAMILTONIAN ENERGETICS FOR 
THE DBL PHASE 

In this section we ask what energetics may stabilize 
the DBL or DLBL phases. Specifically, we motivate and 
study the J — K4 Hamiltonian Eq. ((5]) with competing 
boson hopping and ring exchange terms. When J > 
and K4 > 0, this Hamiltonian does not satisfy the Mar- 
shall sign conditions, so one expects the ground state 
wavefunction to take both positive and negative values. 

The Hamiltonian ([5]) is motivated by considering the 
gauge theory description Eq. (jTTJ) of the DBL phase in the 
strong coupling limit of the gauge theory, h 3> K,t\\,t±. 
In this limit, one can perturbatively eliminate the gauge 
field; the resulting boson Hamiltonian contains, among 
other terms, both Hj and H4 with J = At\\t±/h and 
K4, = A' Ktf, I h A - A"i|^//i 3 with positive numerical 
coefficients A, A' , A". Only the signs of the contributions 
to J and K4 are of interest here, since in what follows 
both J and K4 are taken as free parameters of the Hamil- 
tonian Eq. J8J). Coming from the gauge theory, we have 
J > and K4 > in the regime of primary interest when 
K ~ h and tn ^> t±. To define the system, we also need 
to specify the boson density per site p. The hard-core 
boson model has particle-hole symmetry, so it is enough 
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to consider p < 1/2. 

Before proceeding, we note that the above Hamiltonian 
with J > and < was introduced and analyzed in 
Ref. 0- Since there is no sign problem in this case, exten- 
sive quantum Monte Carlo studies of Refs. ll5lflMl7l were 
able to map out the phase diagram of the model. In addi- 
tion to a bosonic superfiuid phase, the numerical studies 
found two ordered phases at half-filling, a bond-ordered 
stripe phase and a staggered charge density wave (CDW) . 
When the Monte Carlo was performed at fixed chemical 
potential,— varying the chemical potential drove first- 
order transitions out of the commensurate solid phases 
into the superfiuid. On the other hand, when the Monte 
Carlo was performed at fixed density^ the superfiuid 
phase was observed as soon as the density dropped be- 
low half-filling, while a phase separation occurred when 
the density was decreased further. 

The case J > and > was not studied in Monte 
Carlo since there is a sign problem. Below, we present a 
rudimentary variational energetics study of this regime. 
The results are summarized in Fig. [15] and lend some sup- 
port that the DLBL phase of Sec. lIII Cl mav be stabilized 
by such competing hopping and ring exchange terms. 

A. Model with ring exchanges only 

Consider first the Hamiltonian with the ring term only. 
In this case, we can change the sign of K4 by divid- 
ing the square lattice into four sublattices and perform- 
ing a transformation b — > —6 on one sublattice. The 
cited quantum Monte Carlo studies did not consider the 
pure ring model, but extrapolating these to J = 0, it 
seems likely that at half-filling the ground state is the 
(it, 7r) CDW. However, we also expect that there is a 
regime away from half-filling that realizes a novel "Exci- 
tonic Bose Liquid" (EBL) phase discovered and studied 
by Paramckanti et al.j^ which we now discuss. 

The EBL was accessed by considering a rotor version 
of the ring Hamiltonian, 

■ffrotor = — X^("r ~ nf 
r 

- 1^1 X] COS (^r - </>r+x + 0r+i+y - 0r+y) , 
r 

where the phase <f> r and the boson density n r are canon- 
ically conjugate. When the cosine in the ring term is 
expanded, i.e., in the phase with no topological defects, 
one finds a quasi-one-dimensional dispersion relation for 
the excitons (density waves), w 2 . ~ sin 2 (fca;/2) sin 2 (/cj / /2). 
This gapless normal Bose state is expected to be stable 
against becoming an insulator when |/C| dominates over 
U. 

The ring only model on the square lattice has spe- 
cial conservation laws: particle numbers on each col- 
umn and each row are individually conserved. It is 
this property that is responsible for the vanishing of lo^ 



along the lines k x — and k y = in the EBL phase. 
Paramckanti et al^ showed that due to the gapless lines, 
EBL is a critical (power-law) 2d quantum phase with con- 
tinuously varying exponents. 

For our energetics study, it is useful to have a good 
trial wavefunction for the EBL phase of hard-core bosons. 
Motivated by the described "spin wave" theory, we write 
the following wavefunction 

*EBL <0) ( r l> r 2:---) « exp[-^U EB L(ri -Tj)] , 

i<j 

(118) 

with 

( r _ r ') = Iy W e ik.(r-r') 

EBU ' V ^ 4|sin(fc a ./2)||siii(* w /2)| 

(119) 

In the spin wave theory, we have W = y/K/\lC\, but more 
generally W is treated as a variational parameter. The 
summation over k excludes the lines where either k x = 
or k y = 0. Such wavefunction can be defined in any sector 
with fixed number of particles on each row and column 
as is appropriate in the study of the ring-only model. We 
find that for moderate W and away from half-filling, this 
wavefunction has a density structure factor of the form, 
(p(— k)/5(k)) ~ |/ca;||fcy| at long wavelengths as expected 
in the EBL. Furthermore, the EBL theory predicts that 
the box correlator defined in Eq. (O decays as a power 
law with a continuously varying exponent, and this is 
what we observe for the wavefunction. In the energetics 
study, the parameter VV is varied to minimize the ring 
energy H±. 

Returning to the model with if 4 > 0, the trial wave- 
function is 

4BL >0) (ri, r a , . . . ) (x (-1)^ cxp[- £ u EBL (r t - r,)] , 

i<3 

(120) 

where Ni is the number of bosons on the sublattice I 
of the four sublattices. One can readily translate the 
preceding results to this wavefunction, which has the sign 
structure appropriate for the pure ring exchange model 
with K 4 > 0. 

Interestingly, the DLBL wavefunction in the extreme 
D-eccentricity limit, which was introduced towards the 
end of Sec. lIII Cl and discussed further in Sec. lIVB~2l also 
has the correct sign structure for the ring Hamiltonian. 
This follows from an observation that this wavefunction 
takes opposite signs for any two boson configurations that 
are connected by H4. It is important here that H 4 con- 
tains only unit square plackets, and one subtlety is that 
such unit ring moves are not ergodic in the sectors men- 
tioned earlier with fixed boson numbers on each row and 
column. The signs of the EBL and extremal DLBL wave- 
functions as written will not agree everywhere. However, 
in the study of the ring model, it is natural to consider 
smaller sectors comprising states that are connected by 
repeated applications of H4, and in such sectors the signs 
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of the two wavcfunctions agree. Thus, one needs to care- 
fully specify what sector restriction is being made when 
working with these wavcfunctions. As far as the K4 ener- 
getics is concerned, we find essentially no difference be- 
tween the wavefunctions obtained by restricting to the 
smaller (ergodic) sectors or the larger sectors which fix 
only the boson numbers on each row and column. The 
presented results are all for the latter choice of sectors; 
the Monte Carlo sampling is performed using arbitrary 
rectangular ring moves, which are ergodic in this sector. 

Figure [14] shows the ring exchange energy per site for 
the optimized EBL and extremal DLBL wavefunctions 
plotted vs boson density. The two trial energies are 
fairly close; the EBL state has somewhat lower energy 
f° r P ^ 1/4, and the difference is largest near p = 1/2. 
Note that the extremal DLBL wavefunction has no vari- 
ational parameters, while the EBL state has one param- 
eter. The trial DLBL energy can be somewhat improved 
by maintaining the same sign structure but taking a vari- 
able power of the determinants; such trial energies (which 
we do not show) approach closer to the optimal EBL en- 
ergies. Figure Q3] suggests that the DLBL wavefunction 
not only has the right sign structure, but also has simi- 
lar short range correlations to the EBL state. Of course, 
additional Jastrow-type factors can be used to further 
improve the DLBL energy, but here and below we em- 
phasize the good sign structure and correlations that are 
already present in the determinantal construction itself. 

The similarity between the EBL and extremal DLBL 
wavcfunctions is also revealed in their density struc- 
ture factors, see also Sec. IIV B 21 In both states away 
from half-filling, for small \k x \ <C \k y \, we can write 
Db(k x ,k y ) = A(k y )\k x \. In the DLBL state, A(k y ) ap- 
pears to be independent of k y . On the other hand, in 
the EBL state, A{k y ) ~ \k v \. Thus, the character of the 
k x — > singularity is the same in both states for any 
non-zero k y , and it is only in the limit k x ,k y — > that 
the two differ. Interestingly, both wavefunctions also ex- 
hibit "2fci?" singularities in the density structure factor. 
Such singularities are not apparent in the naive contin- 
uum "spin wave" EBL theory, but manifest themselves 
in the lattice hard-core boson EBL wavefunction. 

The study of density correlations also reveals that these 
optimized states at half-filling in fact have spontaneous 
(tt, 7r) charge order, while away from half-filling no or- 
der is observed. Thus, even though in the present crude 
energetics work we do not consider directly conventional 
ordered states of bosons such as charge or bond density 
waves, the EBL and DLBL wavcfunctions provide some 
access to the energetics of such ordered states. (Note, 
however, that these trial states likely have some degree 
of EBL-ness, e.g., finite compressibility; this is analogous 
to the well-known observation that when a conventional 
Jastrow wavefunction is driven into a regime of charge or- 
der, it in fact represents a supersolid that also supports 
off-diagonal long range order.) As we have already men- 
tioned, Monte Carlo studie a 16 ' 17 suggest staggered CDW 
at p = 1/2, and our finding of the same tendency in 




-0.35 



0.1 0.2 0.3 0.4 0.5 

Pbos 

FIG. 14: Trial ring energy in the optimized EBL wavefunc- 
tion (filled circles) and extremal DLBL wavefunction (open 
circles). Measurements are performed on a 24 x 24 system; 
the boson number is always a multiple of 24 since the ex- 
tremal DLBL requires that the number of particles is ex- 
actly the same in each row/column. Away from half-filling 
the optimized EBL and DLBL wavefunctions produce uni- 
form liquids, but exactly at half-filling we find a long-range 
(•7r,7r) CDW order in both states, which is emphasized with 
large box symbols. The thin lines are guides to eye, while the 
thick lines are the result of the Maxwell construction for the 
two trial energies that addresses the possibility of the phase 
separation at low densities. From this construction, the uni- 
form liquids appear to be stable in the the density windows 
0.42 < p < 0.50. 



the considered wavcfunctions lends some support to the 
goodness of the variational approach. Conventional bo- 
son orders arc potentially important near other commen- 
surate boson densities but are less important at generic 
p, and we will not consider charge ordered states further. 

Leaving the charge orders aside, an important feature 
is prominent in the energy plots in Fig. [14] for the con- 
sidered uniform liquids. We see regions where e(p) is 
concave, which signals an instability towards phase sepa- 
ration. Using the Maxwell construction for the available 
data, we conclude that the system phase separates for bo- 
son densities p < 0.42. At low densities, the system will 
spontaneously separate into an empty phase and the EBL 
phase at density p = 0.42. In the region near p = 1/2, 
the considered uniform liquids appear to be stable, ex- 
cept perhaps very close to half-filling. 

The tendency in the ring model towards phase sepa- 
ration at low densities was noted by Rousseau et a/.— 
Indeed, if we consider two bosons on an otherwise empty 
lattice, the lowest ring energy is achieved when the two 
reside on the same square placket, since only in this case 
the ring term is nonzero. Thus, two bosons can be bound 
together by the ring energy. By performing exact diago- 
nalization of small clusters up to 4 x 4 with open bound- 
ary conditions and for different boson numbers, we find 
that the energy per particle is the lowest when the density 
is near half-filling; also, larger open clusters achieve lower 
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energies per boson. One can construct eigenstates of the 
ring Hamiltonian at low densities composed of such iso- 
lated clusters, and the discussed tendencies suggest that 
not only two but many bosons tend to clump together in 
the presence of the ring exchanges. 

One more notable feature in the energy plots occurs 
near half-filling. We show only the p < 1/2 part, while 
the p > 1/2 part would be the mirror image due to 
particle-hole symmetry. Therefore, e(p) shows a cusp 
precisely at half-filling - at least in the EBL case and 
for the system sizes studied. Since \i = de/dp, the cusp 
may be interpreted as an evidence for a charge gap at this 
density, which appears to be in line with the proposal of 
the CDW order by the Monte Carlo studies. 

We do not know whether the above observations persist 
in the limit of large systems. Thus, it may also happen 
that there will be no regime of the EBL phase in the spe- 
cific model if the system chooses to phase separate into 
empty space and a half-filled charge insulator. If there 
is a stable EBL regime over some density window near 
half-filling, then the phase separation at lower densities 
is likely to be into empty space and an EBL phase of an 
appropriate density determined from the Maxwell con- 
struction. The existence of an EBL phase near half-filling 
could presumably be checked by careful Monte Carlo sim- 
ulations for the pure ring model. It is also plausible that 
the regime of the stable uniform EBL phase could be 
extended by adding interactions that disfavor phase sep- 
aration, but we have not explored this in the present 
work. 



B. Full J-K4, Hamiltonian 

We now turn to the case with non-zero boson hopping 
J > 0. Consider first small J, and suppose we start in the 
uniform EBL or DLBL regime assuming such exists. As 
we have discussed, the hopping and ring terms frustrate 
one another. At present, we do not know how to include 
the effects of hopping if we begin with the EBL. Thus, 
the EBL wavefunction is defined in sectors with fixed 
boson numbers on each row/column, but what relative 
signs should we take for the different sectors that are 
now mixed by the boson hopping? If, for example, we 
naively extend Eq. (|120[) to all boson configurations, with 
perhaps modified Jastrow factors, the expectation value 
of Hj remains precisely zero. It seems that one needs to 
change the sign structure of the wavefunction in order to 
treat the frustration accurately. 

On the other hand, the DLBL wavefunction provides 
a natural continuation from the extremal case of the ring 
only model, and we find that it is indeed capable of uti- 
lizing both the ring energy and the hopping energy. More 
precisely, for non-zero J the optimized D-eccentricity is 
finite and varies as J is increased, moving towards zero 
D-eccentricity at large J. Within our restricted varia- 
tional energetics study this is the strongest evidence we 
can offer that the ground state of the ring Hamiltonian 
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FIG. 15: Schematic "variational phase diagram" of the J— i\4 
model. Numerics is done for the same sizes as in Fig. 1141 
The phase boundaries are obtained by examining the DBL 
state, allowing for the possibility of phase separation, and 
the superfluid state (the latter appears to be always uniform 
in our study). The optimal values of t±/t« throughout the 
whole DBL phase are such that the fermions have open Fermi 
surfaces, i.e., this is the DLBL phase of Sec. IIII CI 



in this regime is in the DLBL phase. 

Some details of the variational study are summarized 
in the "phase diagram" in Fig. [131 At a given density, 
the trial energy of the DBL state is optimized using the 
ratio ij_/t|| as a variational parameter. Similarly to the 
pure ring model, we also analyze the possibility of phase 
separation, which is done in the spirit of Fig. [T3J Wc 
find a stable uniform DBL phase in the density range 
roughly similar to the ring only model. Physically, we 
expect the phase separation to be suppressed when the 
hopping becomes non-zero, but our study is too limited 
to explore this systematically. 

The phase diagram Fig. [T5] also includes a conven- 
tional superfluid state of bosons, which is expected in the 
regime of large J. We take the Jastrow-type form for the 
trial wavefunction, Eq. (fTo) . and use the pseudopoten- 
tial u(r — r') = W/\r — r'| p with two variational param- 
eters W and p. Allowing more parameters such as an 
independent nearest-neighbor pseudopotential does not 
visibly change the optimized energies. Interestingly, as 
deduced from the energetics, the superfluid-DBL phase 
boundary roughly coincides with the limit of vanishing 
D-eccentricity t± = tn. The resulting (det) 2 wavefunc- 
tion is positive, and, as we discussed at length in Scc. lIVl 
this state has off-diagonal long range order. Of course, 
it does not have the same long-wavelength correlations 
as the superfluid (e.g., it has unusual 2Uf density cor- 
relations), but it appears to have reasonable short-range 
correlations. 

Finally, we also remark that in the present J — K4 
study, the optimal DBL phase when it wins over the su- 
perfluid is in fact the DLBL phase with open Fermi sur- 
faces discussed in Sec. IIII Cl and is in the regime where the 
boson Green's function decays exponentially in all direc- 
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tions. As we discussed in Sees. Mil and IIV1 despite such 
"locality" in the boson correlator, this unusual phase is 
still critical with many low energy excitations. 

Summarizing, the DBL wavefunction is able to inter- 
polate between the EBL and the superfluid regimes in 
the frustrated J — K$ model. This is the main argument 
supporting our proposal that the ground state of the ring 
Hamiltonian in the intermediate regime is in the DLBL 
phase. 



VIII. SUMMARY AND DISCUSSION 

Our main results were summarized in the Introduc- 
tion, so here we only highlight the most interesting points 
and discuss possible future directions. We are searching 
for examples of uncondensed quantum boson liquids that 
respect time reversal and occur at generic continuously 
varying densities. In this paper, we proposed several 
wavefunctions that produce such liquids and studied in 
detail the specific states dubbed DBL and DLBL. While 
the initial wavefunction construction appears somewhat 
ad-hoc, we argued that it is no more so than the slave 
particle construction of spin liquid states in frustrated 
antifcrromagnets, which has been developed to maturity 
over the last two decades. We introduced a particular 
slave fermion treatment of hard core bosons on a square 
lattice that corresponds to the proposed DBL wavefunc- 
tions. This approach allows to go beyond the trial ground 
state wavefunction, suggesting natural excitations, and 
eventually leads to a gauge theory description of the low- 
energy properties of the putative boson liquid phase. Us- 
ing techniques previously developed for the so-called uni- 
form RVB spin liquid with a Fermi sea of spinons, we 
argued for the possible stability of the new DBL phases 
to gauge fluctuations and found the DLBL regime to be 
particularly stable. Within the gauge theory formulation, 
we also analyzed boson and boson density correlators at 
long distances, which reveal special singular surfaces in 
the momentum space. 

The singular surfaces were confirmed by measuring the 
properties of the proposed DBL /DLBL wavefunctions di- 
rectly; in the same numerics, no potential orders were 
observed, indicating stability of the states. We remark 
here that the Gutzwiller wavefunctions and the gauge 
theory should not be viewed as identical descriptions, 
and, in particular, the detailed long-distance properties 
may be different. The Gutzwiller wavefunctions can be 
used for crude energetics and qualitative characteriza- 
tion, but in our opinion the gauge theory is more com- 
plete as far as the actual quantum phase is concerned. 
Indeed, the wavefunctions as constructed do not incorpo- 
rate any dynamics of the gauge fluxes, which are clearly 
additional low-energy "variational" degrees of freedoms 
that the quantum system is going to explore and utilize. 

From the singular momentum surfaces in the boson 
properties, one can in principle recover the underlying 
Fermi surfaces of the slave particle construction; the lat- 



ter have correct volumes for the boson density, offering 
a tantalizing possibility of some "Luttinger theorem" for 
uncondensed liquids of interacting bosons. 

Some striking thermodynamic properties of the pro- 
posed boson liquids follow from their fermion - gauge 
character ^^M1^^7^9jo In both thc DB L and 

DLBL phases, because of thc gaplcss Fermi surfaces and 
the gauge interactions, the specific heat at low tempera- 
tures is expected to behave as C ~ T 2 / 3 , while the resis- 
tivity is expected to vanish as Rb ~ T 4 / 3 . The last result 
in the DBL phase with closed Fermi surfaces is simply 
a quote from earlier studies of such Fermi sea - gauge 
systemS ) 34 i 35 while the case with open Fermi surfaces re- 
quires slightly more care. Our argument is based on the 
application of thc Ioffc-Larkin rule^ which gives the boson 
resistivity as Rb = Rd t +Rd 2 ■ For the Fermi surfaces as in 
the right panel of Fig. [T] and the conductivity measured, 
say, in the x direction, the small-angle scattering by the 
gauge field can effectively degrade the di fermion current, 
therefore the standard result Rd 2 ~ T 4 / 3 . On the other 
hand, such scattering cannot completely degrade the d\ 
fermion current, so Rd 1 <C Rd 2 , an d Rb ~ T 4 / 3 follows. 
(Parenthetically, the d\ current can be degraded by quar- 
tic or more particle interactions, but such contributions 
vanish as T 2 or faster.) 

Intriguingly, in the DLBL phase, the boson correlators 
are expected to be short-ranged despite such manifest 
thermodynamic signatures of gaplessness and criticality. 
On the other hand, the boson box correlator in the DLBL 
decays as a power law and is negative, — x~ 8 , offering per- 
haps a more physical glimpse of thc fcrmionic partons in 
this phase. Indeed, consider inserting bosons at (0, 0) 
and (x, x) and removing at [x, 0) and (0,x). The domi- 
nant contribution to the box correlator comes from the d\ 
fermions propagating (0,0) — » (x, 0) and (x, x) — » (0,x), 
while the di paths are (0, 0) — ► (0, x) and (x, x) — > (x, 0), 
which is illustrated in Fig. 1161 The minus sign of the box 
correlator is then due to one fermionic exchange needed 
so that the two pairs of injected fermions are removed 
with the same pairing. 

The proposed DBL/DLBL wavefunctions do not sat- 
isfy Marshall signs, and their interesting sign structure is 
brought out by the nodal pictures in the continuum such 
as Fig. O which give some caricature of the wavefunction 
signs also on the lattice even though the nodes are not 
sharply defined in this case. Clearly, these wavefunctions 
cannot be ground states of boson models without frustra- 
tion. To address the question of what Hamiltonians may 
stabilize such phases, we considered a particular frus- 
trated hard core boson model with competing hopping 
and four-site ring exchange terms, and our energetics 
study suggests that the DLBL state is a good candidate 
in this model near half-filling. Unfortunately, the frus- 
trated nature of the boson motion makes this model not 
suitable for large system quantum Monte Carlo studies. 
We still suggest it as an interesting model for numerical 
studies such as exact diagonalization of small systems or 
DMRG; the gapless nature of the boson liquid may make 
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FIG. 16: Boson box correlator (6 1 " (0, 0)6 t (x, x)b(x, 0)6(0, a;)) 
in the DLBL phase reveals the fermionic character of the par- 
tons, decaying as power law and having negative sign from 
the fermion exchange. 



results hard to interpret, but the short-range character of 
the correlations and particularly the quasi-local nature of 
the DLBL phase can perhaps facilitate DMRG to access 
larger systems. 

Looking ahead, we now describe one of the main drives 
behind the study of uncondensed boson liquids, despite 
its own intrinsic interest. We are searching for elec- 
tronic conducting non-Fermi liquids, and the ideas of the 
present work suggest some avenues in this direction. As a 
specific example, consider the slave boson approach that 
is popular in the context of the t — J model of high- 
T c superconductors. The electron operator is written as 
cl = w fl, leading to a theory of spinons and slave bosons 
strongly coupled via an emergent fluctuating gauge field. 
Phcnomenologically, a Fermi sea of spinons is very ap- 
pealing in the strange metal, but what about the slave 
bosons? If they condense, we would recover the conven- 
tional Fermi liquid. It has been argued that the fluc- 
tuating gauge field frustrates the motion of the bosons 



and suppresses the tendency to condense i 4 ' 5 ' 6 ' 34 ' 35 If this 
suppression could persist to zero temperature, we would 
obtain a conducting non-Fermi liquid state. Such a pos- 
sibility is precisely what we are trying to establish in the 
present work, and our approach would be to write the 
slave boson in terms of new (second generation) slave 
fermions, = d\d\. 

On the level of wavefunctions, we would then write an 
electronic state of the form: 

*eiectron(t,l) = [(det) a x (det) v ](t,l) det(f) det(j) , 

where schematically (j) or (J) denotes locations of all 
electrons with spin-up or spin-down respectively. Each 
determinant in the slave boson wavefunction (det) x x 
(det) y evaluates appropriate d-particle orbitals at the lo- 
cations of all electrons irrespective of their spin, which 
assures the no-double occupancy constraint. There is sig- 
nificant freedom in specifying these orbitals, and whether 
such wavefunctions are useful for any electronic system 
requires detailed energetics studies of specific Hamiltoni- 
ans. One can nevertheless develop a low-energy theory 
of such a T = phase in the spirit of the present paper, 
thereby obtaining an itinerant non-Fermi liquid conduct- 
ing phase. If the DBL or DLBL bosonic states studied 
in this work appear as useful caricatures of the behavior 
of the electronic charge, we may even call the resulting 
phase a "D-wave metal"! Such all-fermion description 
of electrons on the square lattice with no double occu- 
pancy has some phenomenological appeal, but detailed 
explorations are left for future work. 
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